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ABSTRACT 


A  time -dependent,  three-dimensional,  large-scale  cloud  model  has  been 
developed  for  the  prediction  of  cloud  cover,  cloud  liquid/ice  water  content 
(LWC/IWC) ,  precipitation,  specific  humidity  and  temperature.  Partial  cloudiness 
is  allowed  to  form  when  large-scale  relative  humidity  is  less  than  100%.  Both 
liquid  and  ice  phases  are  included  in  the  model.  The  liquid  phase  processes 
consist  of  evaporation,  condensation,  autoconversion  and  precipitation.  The  ice 
phase  processes  include  heterogeneous  nucleation  to  generate  ice  crystals, 
depositional  growth  to  simulate  Bergeron-Findeison' s  process,  sublimation  to 
deplete  ice  crystals,  and  gravitational  settling  of  ice  crystals.  The  radiative 
transfer  parameterization  program  is  based  on  a  broadband  method  and  involves  the 
transfer  of  infrared  and  solar  radiation  in  clear  and  cloudy  regions.  The 
broadband  infrared  emissivity,  reflectivity,  and  transmissivity  for  cirrus 
clouds,  as  well  as  the  broadband  solar  absorption,  reflection,  and  transmission 
values  for  low,  middle  and  high  clouds  are  computed  based  on  the  cloud  LWC  and 
IWC  interactively  generated  by  the  cloud  model.  Large  amounts  of  satellite  data, 
including  cloud  cover  climatology  derived  from  the  U.S.  Air  Force  three- 
dimensional  nephanalysis  (3DNEPH)  and  earth  radiation  budget  (ERB) ,  have  been 
processed  into  formats  appropriate  for  verification.  The  96 -hour  model 
simulations  of  cloud  cover,  outgoing  longwave  radiation  (OLR)  ,  and  cloud  LWC  have 
been  verified  against  data  analyzed  from  3DNEPH,  ERB  and  the  Nimbus -7  Scanning 
Multichannel  Microwave  Radiometer  (SMMR)  satellite  observations,  respectively. 
The  predicted  cloud  IWC  is  compared  to  in  situ  observations  as  well  as  to  results 
from  other  studies.  The  predicted  cloud  and  radiation  results  compare  well  with 
those  analyzed  from  satellite  data.  Numerical  experiments  are  carried  out  with 
and  without  radiative  heating  and  ice  phase  processes  in  the  cloud  formation 


scheme.  The  inclusion  of  radiative  heating  produces  a  significant  change  in 
temperature,  cloud  cover  and  total  cloud  water  content,  while  the  inclusion  of 
ice  phase  processes  generates  a  substantial  change  only  in  total  cloud  water 
content.  If  the  cloud  LWC  is  sufficiently  large  to  initiate  effective  Bergeron- 
Findeison's  processes,  the  total  cloud  water  content  decreases,  indicating  that 
gravitational  settling  is  an  efficient  mechanism  in  reducing  the  cloud  IWC. 
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1 .  INTRODUCTION 

Clouds  regularly  occupy  about  50%  of  the  sky.  They  are  the  most  important 
regulators  of  the  earth's  radiation  budget,  i.e.,  the  balance  between  incoming 
solar  radiation  and  outgoing  longwave  radiation  at  the  top  of  the  atmosphere 
(TOA) ,  by  which  the  earth's  climate  is  determined.  On  the  one  hand,  clouds 
reflect  a  large  portion  of  incoming  solar  radiation,  producing  a  cooling  effect 
referred  to  as  the  solar  albedo  effect.  On  the  other,  they  trap  the  outgoing 
longwave  radiation  emitted  from  the  earth,  resulting  in  a  warming  effect  referred 
to  as  the  IR  greenhouse  effect.  Since  the  greenhouse  and  albedo  effects  are 
different  in  sign  as  well  as  magnitude,  the  existence  of  clouds  may  have  a 
profound  impact  upon  the  sensitivity  of  climate  to  external  perturbations,  such 
as  the  Increase  of  atmospheric  carbon  dioxide.  In  particular,  cirrus  clouds, 
which  regularly  cover  about  20%  of  the  globe,  are  optically  thin  and  nonblack. 
The  outcome  of  the  greehhouse-versus-albedo  effects  is  intrinsically  modified  by 
the  radiative  effects  of  nonblack  cirrus  (Liou,  1986). 

Recognizing  the  importance  of  clouds  in  weather  and  climate  processes, 
numerous  modeling  efforts  have  been  carried  out  to  understand  the  relative 
significance  of  solar  albedo  and  IR  greenhouse  effects  involving  clouds.  In 
early  models,  the  formation  of  clouds  was  very  crude  and  lacked  physical 
foundation.  An  arbitrary  fractional  amount  of  cloud  was  assigned  to  a  given 
layer  where  large-scale  condensation  occurs.  All  the  excessive  water  vapor  was 
condensed  and  immediately  precipitated  onto  the  ground.  As  a  result,  there  was 
no  cloud  water  or  cloud  ice  remaining  in  the  clouds.  These  two  elements  are 
necessary  for  the  calculation  of  cloud  radiative  effects.  In  climatic 
perturbation  studies,  Wetherald  and  Manabe  (1988)  assigned  a  cloud  amount  of  80% 
to  condensation  and  found  that  cloud  feedback  processes  enhance  the  sensitivity 


of  the  model  climate  to  an  increase  in  atmospheric  carbon  dioxide.  This 
conclusion  is  in  agreement  with  that  drawn  by  Hansen  et  al.  (1983)  .  The  features 
of  cloud  cover  change  obtained  in  these  studies,  as  well  as  in  that  of  Washington 
and  Meehl  (1984),  are  qualitatively  similar. 

Efforts  have  also  been  focused  on  improving  the  specification  of  the 
radiative  properties  of  clouds  and  cloud  parameterization  schemes  in  GCMs  to 
demonstrate  the  importance  of  cloud- radiation  interaction  to  the  simulation  of 
atmospheric  circulation.  Ramanathan  et  al.  (1983)  demonstrated  that  refinement 
in  radiative  properties  of  non-black  high  cirrus  emissivity  significantly 
improves  the  model  simulation  of  atmospheric  temperatures  as  well  as  that  of 
zonal  mean  winds.  Slingo  and  Slingo  (1988)  reached  the  same  conclusion  in  their 
study  of  cloud  longwave  radiative  forcing.  Using  a  GCM  designed  for  medium  range 
weather  prediction,  Liou  and  Zheng  (1984)  illustrated  the  importance  of  cloud  and 
radiative  processes  to  the  maintenance  of  the  Hadley  circulation. 

In  order  to  incorporate  interactive  cloud  and  radiation  processes  in 
models,  cloud  liquid  water  content  (LWC)  and  ice  water  content  (IWC)  are  required 
for  radiation  calculations.  The  models  must  be  able  to  generate  LWC/IWC.  This 
is  the  so-called  prognostic  method,  which  requires  the  additional  prognostic 
equation(s)  for  the  prediction  of  LWC/IWC.  Sundqvist  (1978)  made  the  first 
attempt  to  develop  a  physically  based  prognostic  large-scale  cloud  scheme  for  use 
in  a  GCM.  In  his  approach,  the  cloud  LWC  is  a  prognostic  variable  determined  by 
various  source  and  sink  terms  which  are  parameterized  in  terms  of  associated  bulk 
quantities.  Also  partial  cloudiness  is  allowed  to  exist  when  the  large-scale 
relative  humidity  is  less  than  100%.  Sundqvist  (1981)  incorporated  the  proposed 
cloud  scheme  into  the  European  Center  for  Medium  Range  Weather  Forecasts  (ECMWF) 
GCM  and  carried  out  a  5-day  integration.  The  forecast  and  cloud  amounts  were 
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compared  with  satellite  photographs,  which  showed  that  the  cloud  patterns  are 
reasonably  well -simulated;  the  predicted  cloud  LUC  is  within  the  observational 
range.  However,  these  studies  did  not  include  radiative  effects  and  ice  phase 
processes.  The  importance  of  including  a  more  realistic  and  physically  based 
parameterization  for  cloud  parameters  in  GCMs  has  recently  been  articulated  by 
Sundqvist  (1988). 

Mitchell  et  al.  (1989)  conducted  a  study  of  cloud  feedback  in  the  UK 
Meteorological  Office  GCM.  A  cloud  liquid  water  budget  equation  was  included  in 
the  model.  However,  ice  clouds  were  not  generated  from  the  budget  equation. 
Water  cloud  radiative  properties  are  interactive  in  the  model  based  on  the 
schemes  developed  by  Liou  and  Wittman  (1979)  for  solar  radiation  and  Stephens 
(1978)  for  IR  radiation.  Roeckner  (1988)  also  included  a  water  budget  equation 
in  the  University  of  Hamburg  GCM  to  undertake  a  study  of  the  feedback  processes 
involving  clouds  due  to  external  radiative  forcing.  Smith  (1990)  incorporated 
a  prognostic  cloud  water  budget  equation  in  a  GCM  and  compared  the  simulated 
cloud  and  radiation  fields  with  satellite  observations. 

Interactive  cloud  formation  in  a  model  setting  would  require  a  detailed 
consideration  of  the  sources  and  sinks  of  water  and  ice  particles.  The  inclusion 
of  separate  empirical  equations  for  cloud  LUC  (e.g.,  Xu  and  Krueger,  1991),  IWC 
(e.g.,  Heymsfield  and  Donner,  1990)  and  precipitation  in  the  GCM  diagnostic 
approach  would  appear  to  be  extremely  difficult  and  inconsistent  with  respect  to 
the  model  physics.  In  order  to  understand  the  role  of  radiative  heating  and  ice 
phase  processes  in  the  formation  of  large-scale  clouds,  we  have  developed  a 
global  stratiform  cloud  model  that  includes  prognostic  cloud  LUC,  IWC  and 
precipitation  equations  as  well  as  an  interactive  radiative  transfer  program. 
The  wind  fields  required  for  the  simulation  are  determined  from  a  GCM.  In 
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Section  2,  we  describe  a  large-scale  cloud  model  and  an  interactive  radiation 
program.  Section  3  contains  bulk  parameterizations  of  various  sources  and  sinks 
for  the  cloud  microphysical  processes  simulated  in  the  model.  Section  4  presents 
comparisons  of  model  predictions  and  verification  data,  including  a  description 
of  model  input  and  verification  data.  Results  on  the  numerical  experiments  on 
atmospheric  radiation  and  ice  phase  processes  are  shown  in  Section  5.  Finally, 
conclusions  are  given  in  Section  6. 

2.  MODEL  DESCRIPTION 

a.  Governing  Eqtiaclons 

We  have  developed  a  time - dependent ,  three-dimensional  cloud-moisture  model 
for  the  prediction  of  water  vapor,  cloud  LWC,  cloud  IWC,  precipitation  and 
temperature.  The  first  four  variables  are  governed  by  the  law  of  mass 
conservation.  The  latent  heat  exchanges  during  the  phase  transitions  of  water 
substances  are  determined  by  the  first  law  of  thermodynamics  through  the 
temperature  field.  The  governing  equations  may  be  written  as  follows: 

*  ^(?T^cos<J»]  .  ^(p^w)  .  pr,w 
=  P(^  *  >  .^Sbf  -  ^  f;  ^  f;  ,  (2.1) 

"  .^(?q,vcos4»)]  . 

-  -  pi,  *  S,  >  F^^  ^  F^'^  ,  (2.2) 

.  .^p^  -  pp  -  sbf  *  F„‘’‘  ♦  f;‘  . 


(2.3) 
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a 


(pqj) 


*  ^(pq,vcos4»]  ^^[pq^Cw-V^)] 


^pQ,  ^  Sb,  -  S^  ^  F*’’’  ♦  F;’ 


(2.4) 


d 


(P)  =  PP 


(2.5) 


In  these  equations ,  an  overbar  denotes  the  average  of  a  variable  over  a  grid 
space.  The  variables  u,  v,  and  w  are  wind  fields  and  are  obtained  from  a  GCM; 

1  and  ^  are  the  longitude  and  latitude,  respectively;  p  is  the  air  density,  T  the 
temperature,  Cp  the  specific  heat  at  constant  pressure,  the  adiabatic  lapse 
rate,  P  the  precipitation  flux,  q^,  qp  q,-  the  mixing  ratios  for  water  vapor, 
cloud  water,  and  cloud  ice,  respectively,  and  L^.,  Ljj  and  the  latent  heat  of 
condensation,  sublimation,  and  melting,  respectively.  Also  in  these  equations, 
Qg,  Qj  and  represent  the  condensational  heating  rates  for  cloud  water,  cloud 
ice,  and  the  radiative  heating  rate,  respectively,  P  is  the  rate  of 
autoconversion  from  cloud  droplets  to  raindrops,  Ej.  is  the  evaporation  of  rain 
water,  and  are  the  rates  of  change  of  the  mass  of  ice  particles  due  to 
Bergeron- Finde ison ' s  process  and  sublimation,  respectively,  and  w^  and  V.|,  are  the 
bulk  terminal  velocities  for  cloud  water  and  ice  particles,  respectively.  In 
addition  the  vertical  eddy  flux  and  horizontal  diffusion  terms  have  been  denoted 
by  Fy  and  FjJ,  respectively,  where  x  can  be  T,  q^,  q^  or  q. . 

The  cloud  cover,  tj  ,  is  derived  based  on  an  average  of  the  specific 
humidities  over  clear  and  cloudy  regions  in  the  form 

q  »  qq^  +  (1  -  n)qo  ,  (2.6) 

where  q^  and  q^  are,  respectively,  the  specific  humidity  in  the  cloudy  and  clear 
regions.  After  rearranging  terms,  the  cloud  cover  is  given  by 
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„  qy/qs  -  qo/qs 
q  “  — = - =z-  • 

qc/qs  -  qo/qs  (2.7) 

Here,  we  have  introduced  the  grid-averaged  saturation  specific  humidity,  q^(T)  , 
which  is  a  function  of  the  grid-averaged  temperature,  T.  For  a  stratiform  cloud, 
the  temperature  within  the  cloud  is  approximately  equal  to  the  grid- averaged 

temperature,  i.e.,  q^  -  q^CT) .  Thus, 

(h  -  ho) 

^  -  (1  -  ho)  ’  (2.8) 
where  h  -  %/%  the  averaged  relative  humidity  and  hg  -  qo/^s  referred  to 
as  the  threshold  relative  humidity,  which  is  a  parameter  to  be  prescribed  to 
close  the  model  equations . 

The  cloud  cover  expression  in  Eq.  (2.8)  denotes  a  linear  relation  between 
the  relative  humidity  and  cloud  cover.  A  similar  empirical  relation  based  on 
observations  has  been  developed  by  Smagorinsky  (1960)  for  use  in  GCMs  (Fig.  1). 
For  the  present  large-scale  studies  we  have  used  the  threshold  relative  humidity 
presented  in  Fig.  1  for  the  calculation  of  high,  middle  and  low  cloud  covers. 

With  the  introduction  of  the  threshold  relative  humidity  (h^  <  1)  ,  partial 
cloudiness  in  the  grid  box  is  allowed  when  the  large-scale  relative  humidity  is 
less  than  100%.  The  cloudy  area  and  clear  region  in  a  grid  box  are  denoted  by 
q  and  (1  -  q),  respectively.  Condensation  occurs  in  the  cloudy  area;  evaporation 
takes  place  in  the  clear  region.  Raindrops  formed  in  a  given  layer  are  not 
subject  to  evaporation  within  that  layer.  Instead,  raindrops  that  fall  through 
this  layer  from  higher  clouds  may  evaporate  in  the  clear  region  of  the  layer. 

Since  the  averaged  terminal  velocity  of  cloud  droplets  is  on  the  order  of 
1  cm/s,  which  is  close  to  the  large  scale  vertical  velocity,  we  would  expect  that 
the  terminal  velocity  of  cloud  droplets  and  vertical  air  motion  roughly  cancel 
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each  other  in  stratiform  clouds.  Therefore,  vertical  advection  of  cloud  LWC  in 
Eq.  (2.3)  is  omitted  in  the  model  calculations.  The  vertical  eddy  fluxes  are 
calculated  from  a  K- theory  approach  (Liou  and  Ou,  1983)  and  the  diffusion  terms 
in  spherical  coordinates  are  parameterized  based  on  the  methodology  described  by 
Washington  and  Williamson  (1977) . 
b.  Cloud  Microphysical  Processes 

The  cloud  microphysical  processes  simulated  in  the  present  cloud  model  are 
illustrated  in  Fig.  2.  Both  the  liquid  and  ice  phase  processes  are  included, 
with  the  liquid  phase  on  the  left-hand  side  and  ice  phase  on  the  right-hand  side 
of  the  diagram.  The  liquid  phase  processes  consist  of  condensation  of  water 
vapor  into  cloud  liquid  water  at  temperatures  warmer  than  -40“ C,  conversion  of 
cloud  liquid  water  into  rain  water  by  means  of  autoconversion,  and  precipitation 
of  rain  water  onto  the  ground.  Evaporation  of  cloud  and  rain  water  Increases  the 
specific  humidity. 

The  ice  phase  processes  include  homogeneous  nucleation  of  water  condensate 
to  cloud  ice  at  temperatures  colder  than  -40"C,  and  conversion  of  cloud  ice  to 
precipitation  ice  by  virtue  of  gravitational  settling,  in  which  the  downward  ice 
flux  is  calculated  by  integrating  the  mass  and  fall  velocity  of  an  individual  ice 
crystal  over  the  whole  spectrum  of  ice  crystals.  The  sublimation  of  cloud  ice 
increases  the  specific  humidity.  At  temperatures  above  the  freezing  point,  the 
melting  process  converts  cloud  ice  to  cloud  water. 

An  important  link  between  the  liquid  and  ice  phases  is  the  so-called 
Bergeron-Findeison  process.  The  saturation  vapor  pressure  over  ice  is  less  than 
that  over  water.  As  a  result,  ice  crystals  grow  by  diffusion  at  the  expense  of 
the  supercooled  liquid  cloud  droplets.  Once  generated  by  ice  nucleation,  ice 
crystals  grow  by  deposition  due  to  Bergeron-Findeison' s  process  and  are  depleted 
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by  sublimation  as  well  as  the  gravitational  settling.  Above  the  freezing  point, 
clouds  consist  entirely  of  cloud  water.  In  the  region  0°  to  -40*C,  cloud  water 
and  cloud  ice  may  coexist.  At  temperatures  below  -40° C,  cloud  ice  dominates 
because  of  homogeneous  nucleation. 

The  formation  of  rimed  crystals,  such  as  graupel  and  hail,  is  not 
considered  in  this  model  since  they  occur  more  frequently  in  convective  clouds 
than  in  large-scale  stratiform  clouds, 
c.  Numerical  Implementation 

In  the  horizontal,  the  grid  point  space  is  characterized  by  48  equally- 
spaced  longitudinal  points  and  38  Gaussian  latitudes,  which  is  the  same  as  the 
horizontal  grid  structure  used  in  the  Air  Force  Geophysics  Laboratory  (AFGL)  GCM 
(Yang  et  al .  ,  1989).  We  have  used  this  model  to  provide  the  present  cloud  model 
with  initial  data  and  a  wind  field.  The  horizontal  grid  space  in  the  model  has 
been  modified  near  the  poles  to  maintain  nearly  constant  geographical  distances 
oetween  the  grid  points.  In  the  vertical,  a  16-layer  stretched  z-coordinate  is 
used  in  the  cloud  model. 

The  upstream  scheme  is  utilized  as  a  means  of  numerically  differencing  the 
horizontal  and  vertical  advection  terms.  The  scheme  is  conditionally  stable  and 
its  stability  criterion  follows  the  Courant-Friedrichs-Levy  (CFL)  condition, 
i-e.,  <  1,  where  is  the  maximum  velocity  in  the  computational 

domain.  At  is  the  time  step  and  AS  is  the  grid  size.  This  criterion  is  easily 
satisfied  in  large  scale  models.  For  example,  in  the  vertical  where  w  —  10® 
cm/s.  At  -  10^,  and  AS  -  10^  cm,  the  criterion  is  much  less  than  1.  The  time 
step  used  in  this  study  is  30  minutes.  The  radiation  calculations  are  updated 
every  3  hours,  i.e. ,  the  radiative  heating  rates  are  calculated  every  3  hours  at 
every  horizontal  grid  point. 
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d.  Radiative  Transfer  Parameterization 

The  radiative  transfer  parameterization  program  used  in  this  study  is  based 
on  a  broadband  method  and  involves  the  transfer  of  thermal  IR  and  solar  radiation 
in  clear  and  cloudy  regions.  In  a  clear  atmosphere,  the  entire  IR  spectrum  is 
divided  into  five  bands;  three  for  HjO,  one  for  COj,  and  one  for  O3  absorption. 
The  parameterizations  of  these  broadband  IR  emissivities  were  developed  by  Liou 
and  Ou  (1981)  and  Ou  and  Liou  (1983).  The  solar  spectrum  consists  of  25  bands: 
six  for  HjO,  one  for  CO2,  which  overlaps  the  HjO  2.7-pm  band,  and  18  for  O3. 

In  a  cloudy  atmosphere,  low  and  middle  clouds  are  treated  as  blackbodies 
in  the  IR  radiative  transfer  calculation.  The  broadband  IR  emissivity, 
reflectivity,  and  transmissivity  for  high  clouds,  as  well  as  the  broadband  solar 
absorption,  reflection,  and  transmission  values  for  low,  middle,  and  high  clouds , 
were  computed  based  on  the  cloud  LWC  and  IWC  interactively  generated  by  the  cloud 
model.  The  cloud  radiative  properties  were  calculated  based  on  the 
parameterizations  developed  by  Liou  and  Wittman  (1979).  Accuracy  of  all  the 
aforementioned  parameterization  was  verified  via  more  comprehensive  and  exact 
radiative  transfer  calculations  described  by  Ou  and  Liou  (1988) . 

The  radiative  heating  rate  at  level  z  is  related  to  the  divergence  of  the 
net  fluxes  and  is  given  by 


.  -J_(AF/Az) 

^  pCp 


(2.9) 


where  Uz  is  the  model  layer  thickness  and  AF  is  the  net  radiative  flux  difference 
between  the  layer  top  and  bottom.  This  equation  applies  to  both  IR  and  solar 
heating  rate  calculations  with  downward  solar  flux  and  upward  IR  flux  defined  as 
positive.  For  the  transfer  of  solar  radiation  in  a  cloud  layer,  the  net  flux 
density  decreases  from  the  cloud  top  to  the  bottom  because  of  cloud  absorption 
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and  scattering.  The  absorbed  radiant  energy  leads  to  the  heating  of  the  cloud 
layer.  For  black  clouds  (low  or  middle  clouds),  the  IR  radiative  flux  is 
proportional  to  the  fourth  power  of  the  absolute  temperature  of  the  cloud.  Since 
the  atmospheric  temperature  usually  decreases  with  height,  the  downward  IR  flux 
emitted  by  the  atmosphere  above  a  black  cloud  top  is  much  smaller  than  that 
emitted  upward  by  the  black  cloud.  As  a  result,  the  net  flux  near  the  cloud  top 
is  greater  than  zero  and  leads  to  strong  IR  cooling  in  the  region.  Similarly, 
there  is  IR  heating  near  the  cloud  bottom. 

Specification  of  the  overlap  that  occurs  between  cloud  layers  is  required 
for  the  calculation  of  radiative  transfer  in  a  radiation  model.  In  this  study, 
the  model -generated  clouds  are  strapped  into  three  cloud  decks.  Low,  middle,  and 
high  clouds  fill,  respectively,  layers  (4,5,6),  layers  (7,8),  and  layers  (9,10). 
The  cloud  cover  for  each  deck  of  clouds  is  obtained  by  averaging  model  cloud 
covers  as  follows; 

lie  =  (n4  +  »l5  116)  /  3  ,  (2.10a) 

=  (»J7  +  tia)  /  2  ,  (2.10b) 

%  °  (^9  ’’lio)  /  2  .  (2.10c) 

where  and  denote  deck  cloud  covers  for  low,  middle,  and  high  clouds, 

respectively. 

To  determine  the  total  cloud  amount,  we  have  assumed  that  clouds  overlap 
each  other  in  a  statistically  random  manner.  Thus  the  total  cloud  cover  is  given 

q  =  1  -  (1  -  qt)(l  -  t^„)(l  -  q^)  .  (2.11) 

where  q  is  the  total  cloud  cover  over  a  grid  area.  For  partly  cloudy  conditions, 
the  radiative  heating/cooling  rate  at  each  model  layer  is  obtained  by  linearly 
weighting  the  percentages  of  the  total  cloud  cover,  q,  and  clear  portion,  (1  - 
q ) ,  in  the  form 
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(11)1*  =  TI(|I)=  -K  (1  -  .  (2.12) 

where  the  superscripts  pc,  c,  and  nc  represent  partly  cloudy,  cloudy,  and  clear 
conditions,  respectively. 


3.  BULK  PARAMETERIZATION  OF  CLOUD  MICROPHYSICS 


Since  it  is  not  practical  in  a  large-scale  cloud  model  to  simulate  detailed 
microphysical  processes  with  respect  to  each  individual  cloud  particle,  we  must 
relate  these  processes  in  terms  of  bulk  quantities,  which  represent  the  various 
sources  and  sinks  for  cloud  water  substances, 
a .  Liquid  Phase 
(1)  Condensation 

Condensation  occurs  as  a  result  of  expanding  air  under  saturation 
conditions  and  may  be  related  to  the  time  rate  of  change  of  the  saturation  vapor 
pressure.  Consequently,  the  condensational  heating  rate  may  be  derived  from  the 
first  law  of  thermodynamics,  the  Claus ius-Clapeyron  equation  and  the  hydrostatic 
approximation  as  follows: 


Qc  • 


qs(T) 

“RT 


LR,  -  CpR^T  ' 

CpR^T^  +  ^l2. 


w  , 


(3.1) 


where  is  the  dry  adiabatic  lapse  rate,  is  the  gas  constant  for  air,  is 
the  gas  constant  for  water  vapor,  L  is  the  latent  heat  for  liquid  and  vapor 
phases,  and  q^  is  the  saturation  vapor  pressure.  The  condensation  is  strongly 
dependent  on  the  vertical  velocity,  w.  As  saturated  air  rises,  i.e.,  w  >  0, 
condensation  takes  place  and  cloud  forms,  while  depletion  of  clouds  occurs  with 
descending  saturated  air. 


(11)  Autoconversion 
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Autoconversion  is  a  process  through  which  small  cloud  droplets  merge  into 
large  raindrops.  The  parameterization  of  autoconversion  follows  the  simple 
exponential  form  developed  by  Sundqvist  (1978); 

P  =  Coqid  -  exp[-(qi/(Tiqir))2]  )  ,  (^-2) 

where  Cq  denotes  the  typical  time  for  the  conversion  of  droplets  to  raindrops, 
and  q^^  is  a  reference  value  of  cloud  LWC  above  which  the  conversion  from  cloud 
droplets  to  raindrops  increases  rapidly. 

(lii)  Evaporation  of  raindrops 

The  evaporation  of  raindrops  may  be  derived  from  the  combination  of  the 

Pick  law  of  mass  diffusion  and  the  Marshall -Palmer  raindrop  size  distribution. 

The  size  distribution  is  assumed  to  remain  the  same  throughout  the  whole  process 
of  evaporation.  The  raindrop  evaporation  rate  is  given  by 

Ep  =  Ke(l  -  h)P°*^2  (3.3) 

where  the  evaporation  rate  constant  is  1.785  x  lO'^q^^. 
b.  Ice  Phase 

The  ice  phase  parameterizations  include  homogeneous/heterogeneous 
nucleation  to  generate  ice  crystals,  depositional  growth  to  simulate  the 
Bergeron-Findeison  process,  sublimation  and  melting  of  ice  crystals,  and 
gravitational  settling  to  deplete  the  ice  crystals.  The  Bergeron-Findeison 
process  occurs  primarily  at  temperatures  between  0*  and  -40*C,  while  below  -40“C 
homogeneous  nucleation  dominates.  The  parameterizations  of  these  source  and  sink 
terms  are  described  below. 

(i)  Homogeneous  nucleation 

Homogeneous  nucleation  takes  place  whenever  condensation  occurs  at 
temperatures  below  -40*C.  Homogeneous  nucleation  is  analogous  to  the  liquid 
water  condensation  process,  and  may  be  written  in  the  form 
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Qi  '  tlTrf 


qs(T) 


LdRa  -  CpRvT 


w  , 


(3.4) 


(  CpR^T2  +  CkL^J 
where  is  the  latent  heat  for  ice  and  vapor  phases. 

(ii)  Heterogeneous  nucleation 

Heterogeneous  nucleation  leads  to  the  formation  of  ice  crystals  on  active 
ice  nuclei  primarily  at  temperatures  between  0*  and  -40®C.  Based  on  the 
experimental  work  of  Mason  (1971)  and  Koenig  and  Murray  (1976),  it  is  assumed 
that  the  active  ice- forming  nuclei  (IN)  concentration  may  be  expressed  in  terms 
of  temperature  as  follows : 

N  =  Aiexp|-Ill^max[(T  -  Tq)  ,  T*]| 


(3.5) 


where  -  1,  Aj  -  4,  T*  <-  O'C,  and  Tq  is  a  threshold  temperature,  which  is  taken 
to  be  -20* C.  Here,  it  is  assumed  that  when  the  saturated  air  temperature  is 
lower  than  Tq,  the  total  number  of  IN  remains  the  same  as  that  at  Tq. 

The  heterogeneous  nucleation  rate  of  ice  crystals,  Sj,,  is  given  by 

-  n,  dN 

^  ^■3t  ’  (3.6) 


where  iDq  is  the  mass  of  a  newly  nucleated  ice  crystal,  prescribed  as  10"'^  grams, 
which  is  equivalent  to  an  ice  sphere  with  a  diameter  of  2.77  pm  (Koenig  and 
Murray,  1976),  and 


dN 

■3t 


if  T  i  -20*C 
St 

0  if  T  <  -20"C  . 


(3.7) 


Here,  N^^^^  and  N^  are  the  IN  concentrations  at  consecutive  time  steps  of  t  +  At 
and  t,  respectively. 

(iii)  Bergeron-Findeison  process 
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The  Bergeron-Findeison  process  transforms  cloud  water  into  cloud  ice  at 
temperatures  between  0°  and  -40°C.  To  simplify  the  complicated  vapor  diffusion 
equation,  Koenig  (1971)  developed  a  simple  equation  to  describe  the  crystal  vapor 
diffusion  growth  rate  as  a  function  of  mass.  The  equation  is  written  as  follows: 


dm  *4 


(3.8) 


where  aj  and  a^  are  temperature -dependent  parameters  given  in  Koenig  (1971) ,  and 
m  is  the  crystal  mass  in  grams. 

The  bulk  quantity  of  the  ice  crystal  depositional  rate,  denoted  as  Sgp ,  may 
be  derived  from  Eq.  (3.8)  by  multiplying  both  sides  by  N,  as  follows: 


Sb,  =  .^(pq,)  =  tif,N’^S(pqi)“^ 


(3.9) 


where  fj  -  9i/(qi  +  q,)  is  the  cloud  liquid  fraction  serving  as  an  adjustment 
factor,  which  allows  for  higher  ice  crystal  growth  rates  with  more  cloud  water. 
The  growth  is  terminated  if  there  is  no  cloud  water  to  supply  water  vapor, 
(iv)  Sublimation  and  melting  of  ice  crystals 

Similar  to  the  derivation  for  the  evaporation  of  raindrops,  we  have 
developed  the  following  equation  for  the  rate  of  sublimation: 


Ss  =  0.1037  (pq^)  (1  -h)  (pqi)°-*5 


(3.10) 


At  temperatures  above  0*C,  cloud  ice  crystals  are  assumed  to  melt 
instantaneously. 

(v)  Gravitational  settling 

Gravitational  settling  is  an  important  sink  for  cloud  IWC.  Therefore,  it 
is  desirable  to  have  a  representative  (mean- volume -weighted)  downward  flux  of  ice 
water.  Following  Starr  and  Cox  (1985),  the  downward  ice  flux  may  be  expressed 
as 
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VxCpqi)  =  nmv  dL  , 


(3.11) 


where  n,  m,  and  v  are  the  ice  crystal  size  distribution,  mass  and  terminal 
velocity  of  an  ice  crystal,  respectively,  and  L  is  the  maximum  dimension  of  an 
ice  crystal.  and  are  the  minimum  and  maximum  sizes  of  ice  crystals. 
The  cloud  IWC  may  be  further  related  to  n  and  m  as 


pqj  =  [‘■"“’‘nmdL  . 


(3.12) 


In  order  to  calculate  the  downward  ice  flux,  the  size  distribution,  mass,  and 
terminal  velocity  for  the  precipitating  ice  crystals  need  to  be  specifically 
expressed  in  terms  of  the  maximtom  dimension  of  ice  crystals. 

The  size  distributions  of  ice  crystals  in  mid- latitude  cirrus  clouds  have 
been  measured  from  aircraft  using  an  optical  probe  (Heymsfield,  1977).  Using 
these  data,  Heymsfield  and  Platt  (1984)  have  averaged  all  spectra  of  crystals 
within  every  5'C  temperature  interval  ranging  from  -20*C  to  -60*C.  The  ice 
crystal  size  distributions  in  terms  of  the  maximum  dimension  of  ice  crystals  may 
then  be  parameterized  for  different  temperature  ranges  in  the  form 

n  =  ,  (3.13) 


where  a^  and  b^,  are  temperature  dependent  empirical  coefficients. 

The  minimum  length  that  the  optical  probe  can  measure  is  about  20  pm. 
Small  ice  crystals  less  than  20  pm  in  size  could  be  missed  by  the  measurement 
technique.  Consequently,  small  ice  crystal  data  were  excluded  from  the  analysis 
by  Hejmsfield  and  Platt  (1984) .  For  this  reason,  20  pm  is  taken  as  the  minimum 
length  of  ice  crystals,  l^j„,  which  is  needed  in  the  calculations  of  Eqs.  (3.11) 
and  (3.12). 
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(vl)  Crystal  mass  and  terminal  velocity 

The  individual  ice  crystal  mass,  m,  may  be  parameterized  in  terms  of  the 
maximum  dimension  of  an  ice  crystal,  L,  as  follows: 


m  = 


(3.14) 


where  m  and  L  are  in  units  of  grams  and  pm,  respectively.  The  constants  a^ji  and 
associated  with  various  ice  crystal  habits  are  given  in  Heymsfield  (1972). 
The  terminal  velocity  v  in  m/s  of  an  ice  crystal  of  maximum  dimension,  L, 
in  pm  is  given  by 

*  (3.15) 

where  the  constants  a^  and  b^  are  dependent  on  the  ice  crystal  habit  and  size  as 
given  in  Starr  and  Cox  (1985) .  These  coefficients  were  derived  at  an  ambient 
pressure  of  400  mb.  For  an  ice  crystal  at  a  given  pressure  p.  Beard  and 
Pruppacher  (1969)  have  suggested  that  its  terminal  velocity  may  be  expressed  by 

v(p)  =  v(p,)(pyp)’/5  ,  (3.16) 

where  p,  denotes  a  reference  pressure  level,  and  v(p*)  is  the  terminal  velocity 
of  a  crystal  at  the  reference  level.  To  obtain  the  terminal  velocity  at  any 
given  pressure,  we  simply  apply  the  pressure  adjustment  to  the  expression  in  Eq. 
(3.15). 

(vii)  Downward  ice  flux 

The  ice  crystal  size  distribution,  ice  crystal  mass,  and  terminal  velocity 
are  all  parameterized  in  terms  of  the  maximum  dimension.  Thus,  the  downward  ice 
flux  may  be  obtained  by  analytically  integrating  Eq.  (3.11),  provided  that  the 
size  of  the  largest  crystal,  is  known.  Since  the  lower  limit,  is  set 

to  be  20  pm,  as  stated  previously,  the  upper  limit  (1^^)  may  be  determined  from 
Eq.  (3.12)  as  follows:  We  first  obtain  pq^  for  the  nth  time  step.  Subsequently, 
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Eqs.  (3.13)  and  (3.14)  are  substituted  into  Eq.  (3.12)  to  calculate  Having 
determined  the  downward  ice  flux  can  then  be  computed  from  Eq.  (3.11)  in 
the  form 

b„ « b. .  b,  ♦  1  ■'-•I"  '■  ‘  ' 

The  downward  flux  is  a  function  of  which  in  turn  is  a  function  of  cloud 
IWC. 

4.  MODEL  PERFOEIMANCE  AND  VERIFICATION 
a.  Model  Input  and  Verification  Datasets 

The  atmospheric  data  generated  from  the  AFGL  spectral  global  model  (Yang 
et  al .  ,  1989)  analyzed  at  12Z  July  1-4,  1979,  were  used  as  inputs  to  the  cloud 
model.  This  period  has  been  chosen  because  of  the  availability  of  cloud  and 
radiation  budget  data  that  have  been  analyzed  on  a  daily  basis.  The  input  data 
contain  temperature,  specific  humidity,  and  wind  fields  (u,  v,  w)  at  all  12  o- 
layers  in  the  atmosphere,  plus  the  surface  pressure  and  geopotential  height. 
These  analyses  were  carried  out  using  a  rhomboidal  15  (R15)  truncation  which 
corresponds  to  a  horizontal  resolution  of  38  and  48  grid  points  in  latitude  and 
longitude,  respectively.  In  the  vertical,  since  the  GCM-generated  input  data  are 
available  only  at  the  o-coordinate ,  it  is  necessary  to  perform  vertical 
interpolations  to  obtain  these  data  in  the  z-coordinate . 

The  3DNEPH  cloud  data  sets  have  been  determined  from  various  satellite, 
aircraft  and  ground-based  observational  sources.  Based  on  3DNEPH  analysis,  the 
global  cloud  climatologies  for  January  and  July  1979  have  been  processed  into 
appropriate  formats  for  model  verification  (Koenig  et  al . ,  1987).  The  global 
distribution  of  cloud  amounts  on  July  1-4,  1979  is  derived  from  the  above  cloud 
climatologies  and  used  as  the  verification  dataset  for  the  cloud  cover  simulated 
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in  the  model . 

It  is  important  to  recognize  the  limitation  of  3DNEPH  cloud  analysis.  A 
threshold  method  is  employed  to  process  satellite  images  to  determine  cloud  cover 
by  a  threshold  temperature.  If  an  observed  brightness  temperature  of  a  pixel  is 
colder  than  the  threshold  temperature,  then  a  cloudy  condition  is  assigned  to  the 
pixel.  Consequently,  the  cloud  amount  tends  to  be  overestimated  when  the 
associated  surface  temperature  is  low,  and  underestimated  when  a  temperature 
inversion  is  located  in  the  region.  This  is  evident  from  the  fact  that  3DNEPH 
over-  and  underestimates  cloud  amounts  in  the  wintertime  and  summertime  polar 
regions,  respectively  (Henderson-Sellers ,  1986).  However,  several 
intercomparisons  of  the  cloud  cover  between  3DNEPH  and  the  other  cloud  retrieval 
techniques  have  been  made  (Koenig  et  al.,  1987;  Hughes,  1984).  The  3DNEPH  total 
cloud  cover  appears  to  be  reliable,  except  in  the  polar  regions  where  all 
satellite  retrieval  methodologies  tend  to  fail.  Additionally,  the  surface 
observational  cloud  climatology  compiled  by  London  (1957)  is  used  as  a 
complementary  verification  dataset  for  the  polar  regions. 

The  archived  ERB  flux  data  tapes  for  July,  1979,  provided  by  the  NASA  ERB 
team  (Kyle  et  al. ,  1990),  are  used  to  validate  OLR  produced  by  the  model.  Since 
ERB  fields  are  strongly  modulated  by  the  radiative  effects  of  clouds,  including 
cloud  cover  and  cloud  LWC/IWC,  they  have  been  widely  used  in  large-scale  cloud 
model  verifications  (Sllngo,  1987;  Smith,  1990).  Monthly  mean  distributions  of 
liquid  water  have  been  derived  from  the  Nimbus  7  SMMR  observations  over  the 
oceans  for  the  period  November  1978  to  November  1979  (Prabhakara  and  Short, 
1984) .  The  liquid  water  estimate  for  July  1979  is  applicable  to  the  monthly 
average  over  an  area  approximately  3*  by  5* .  The  liquid  water  derived  from 
satellite  microwave  channels  includes  vertically  integrated  amounts  of  liquid 
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droplets  in  both  clouds  and  rain.  The  liquid  water  verification  data  are 
confined  to  the  oceans  equatorward  of  60“  latitude.  Over  land  and  ice  surfaces, 
due  to  the  limitation  of  microwave  techniques,  there  is  no  climatological  liquid 
water  data  of  global  extent  for  comparisons. 

Since  there  is  no  ground  truth  with  which  to  compare,  the  accuracy  of 
liquid  water  observations  from  SMMR  is  roughly  estimated  as  100  g/m^,  based  on 
theoretical  considerations  (Prabhakara  et  al.  ,  1983)  .  However,  by  taking  monthly 
averages,  this  accuracy  may  be  improved  to  50  g/m^  (Prabhakara  and  Short,  1984). 
b.  Design  of  the  Verification  and  Initial  Data 

In  the  development  of  any  parameterization  scheme,  the  prediction  results 
must  be  carefully  verified  to  assess  its  performance.  The  large  scale  cloud 
model  was  integrated  for  96  hours  with  a  time  step  of  30  minutes  from  the  initial 
conditions  taken  at  12Z  on  July  1,  1979.  The  wind  fields  are  prescribed  from  a 
GCM.  It  is  noted  that  cloud-radiation  feedback  to  the  dynamic  structure  is  not 
accounted  for  in  this  study. 

With  input  of  the  initial  data  and  wind  fields,  the  cloud  model  is  capable 
of  forecasting  large-scale  cloud  cover  and  cloud  LWC/IWC,  which  are  interactively 
used  in  the  radiation  model.  The  radiation  model  computed  OLR  for  comparison 
with  the  ERB  satellite  data.  In  this  study,  the  interactive  radiation 
computations  are  undertaken  every  3  hours,  while  the  calculated  heating/cooling 
rates  are  input  in  the  cloud  model  at  every  time  step  during  those  3  hours. 

The  simulated  cloud  cover,  OLR,  cloud  LWC  and  cloud  IWC  are  verified 
against  3DNEPH  cloud  cover,  ERB  data,  SMHR  results,  and  in  situ  observations  and 
other  model  studies,  respectively.  However,  there  are  several  problems  which 
make  the  use  of  satellite  observations  in  model  verification  not  entirely 
straightforward.  For  instance,  verification  datasets  such  as  3DNEPH  and  ERB  are 
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based  on  the  daily  average  over  local  time,  while  the  output  of  the  model 
simulation  is  the  instantaneous  global  distribution.  Allowing  intercomparisons 
to  be  made  between  the  satellite  observations  and  model  simulations,  we  have  to 
average  the  model  predicted  fields  over  local  time  in  order  to  obtain  daily 
averaged  model  results. 

The  analysis  data  on  12Z,  July  1,  1979  for  temperature  and  specific 

humidity  are  used  as  initial  values  to  perform  a  96-hour  model  simulation.  Data 
for  the  cloud  LWC/IWC  are  not  available  for  Initial  values  on  a  global  basis. 
At  this  point  we  do  not  have  any  physical  means  to  provide  the  initial  cloud 
LWC/IWC  values.  For  this  reason,  we  have  set  these  values  as  zero,  initially. 
Unless  specified,  all  of  tne  results  including  satellite  observations  and  model 
simulation  are  averaged  over  July  3-4,  1979.  That  is,  most  of  the  results  shown 
in  this  and  latter  sections  are  presented  as  a  2-day  average  of  July  3  and  4. 
Due  to  the  limitation  of  the  initial  cloud  LWC/IWC  fields,  the  results  for  July 
1  and  2  have  not  been  included  in  the  analysis  and  verification, 
c.  Model  Results  and  Verification 

In  Fig.  3a,  the  model-predicted  zonal  mean  total  cloudiness,  averaged  over 
July  3-4,  is  compared  with  the  corresponding  3DNEPH  cloud  data.  Differences 
between  the  two  curves  are  generally  within  10%  of  the  total  cloud  cover  except 
in  the  tropics,  Arctic,  and  Antarctic  areas. 

As  mentioned  in  Subsection  4a,  3DNEPH  under-  and  overestimates  cloud  covers 
in  Arctic  and  Antarctic  regions,  respectively,  during  the  Northern  Hemisphere 
summer.  Henderson-Sellers  (1986)  compared  3DNEPH  monthly  mean  data  with  London's 
(1957)  cloud  climatology  which  is  based  on  surface  cloud  observations.  Since 
London's  cloud  climatology  was  compiled  only  for  the  Northern  Hemisphere,  the 
comparison  is  confined  to  the  Northern  Hemisphere.  London's  cloud  climatology 
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shows  about  70%  of  total  cloudiness  in  the  summertime  Arctic  region  and  40%  in 
the  wintertime.  The  70%  of  total  cloudiness  in  the  Arctic  is  close  to  our  model 
prediction  in  the  area.  This  explains  the  differences  between  the  predicted  and 
3DNEPH  cloud  cover  in  the  Arctic  and  Antarctic  regions. 

The  observed  large  cloudiness  in  the  tropics  is  largely  caused  by  tropical 
anvils,  which  originate  from  convective  clouds.  Since  the  present  large-scale 
cloud  model  has  been  developed  primarily  for  stratiform  clouds,  the  simulated 
cloud  cover  in  the  tropics  is  underestimated,  as  is  evident  in  Fig.  3a.  Another 
shortcoming  in  the  simulated  cloud  cover  is  the  overestimation  of  cloud  cover 
near  60*  S,  where  the  absorbed  solar  radiation  becomes  smaller.  As  a  result, 
strong  radiative  cooling  is  produced  at  the  top  of  middle  clouds,  leading  to  a 
large  increase  of  total  cloudiness  in  the  region.  It  appears  that  an  improvement 
could  be  made  if  a  dynamic  instability  feedback  were  allowed  in  the  simulation. 
The  simulated  small  cloud  cover  loc*>*--d  at  the  subtropic  highs  in  both 
hemispheres  is  realistic.  Ln-^ge  cloudiness  in  the  tropics  is  also  predicted, 
although  the  3DNEPH  cloud  data  suggests  that  the  predicted  cloud  cover  is 
underestimated. 

Figure  3b  reveals  good  agreement  in  OLR  between  the  model  simulation  and 
ERB  observation,  with  differences  generally  less  than  10  W/m^,  except  in  the 
tropics.  The  OLR  in  the  tropics  predicted  from  the  model  is  30  W/m^  highe>  than 
that  analyzed  from  the  ERB  data.  This  overestimation  of  OLR  may  be  explained  by 
the  following  two  factors.  First,  the  altitude  of  high  clouds  defined  in  the 
model  is  usually  lower  than  the  actual  height  of  tropical  anvils.  Second,  the 
simulated  cloud  cover  in  the  tropics  is  underestimated.  Both  curves  in  Fig.  3b 
indicate  two  maxima  located  in  the  subtropics  in  both  hemispheres  with  a  minimum 
in  the  tropics.  These  two  well-defined  OLR  maxima  are  associated  with  the 
subtropical  highs,  where  small  cloudiness  coupled  with  high  temperature  produces 
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large  OLR.  The  minimum  OLR  in  the  tropics  is  related  to  the  ITCZ  region,  where 
thick  high  clouds  are  abundant. 

Data  on  cloud  LWP  are  very  limited.  There  ar3  no  observational  LWP  data 
on  a  daily  basis  with  which  to  compare,  nor  are  the  zcnal  averaged  values  of  LWP 
available  for  July  of  1979.  Since  the  time  domain  in  this  study  is  the  very 
beginning  of  July,  we  have  compared  the  predicted  zonal  mean  LWP  with  the  monthly 
averaged  zonal  mean  data  analyzed  from  SMMR  for  June,  1979.  Figure  3c  shows  the 
zonal  mean  LWP  from  model  simulations  and  from  ^iiMR  measurements  taken  from 
Prabhakara  and  Short  (1984).  The  zonally  averaged  LWPs  over  the  oceans  are  shown 
in  the  figure.  The  zonal  mean  simulated  cloud  LWPs  compare  quite  well  with  SMMR 
observations.  Differences  between  the  two  are  within  50  g/m^  which  is  the 
uncertainty  of  SMMR  data  suggested  by  Prabhakara  and  Short. 

The  maximum  LOT  that  occurs  in  the  tropics  is  related  to  the  ITCZ.  The  two 
well-defined  local  LWP  maxima  at  mid-latitudes  in  both  hemispheres  are  associated 
with  storm  tracks.  The  minimum  LWP  is  found  in  regions  of  subtropical  highs  with 
prevailing  downward  motions.  The  SMMR  data  show  maximum  values  of  about  140,  47, 
and  40  g/m^  in  the  tropics  and  mid- latitudes,  respectively.  The  corresponding 
values  obtained  from  simulation  are  130,  53,  and  65  g/m^,  respectively.  The 
ratio  of  maximum  observed  LWP  in  the  tropics  to  that  in  the  middle  latitude  is 
about  3  to  1.  This  ratio  is  about  2  to  1  based  on  the  model  simulation.  This 
difference  indicates  that  the  model  underestimates  the  cloud  LWP  in  the  tropics, 
where  a  significant  amount  of  LWP  observed  by  the  SMMR  is  associated  with 
convective  clouds,  which  are  not  simulated  in  the  present  model. 

The  predicted  cloud  cover,  OLR  and  LWP  in  the  following  figures  are 
presented  in  terms  of  the  horizontal  cross-sections  to  demonstrate  the 
geographical  distribution  of  cloud  fields.  Cloud  cover  and  OLR  are  illustrated 
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by  the  gray  shades  in  which  the  intensity  in  the  picture  is  proportional  to  a 
certain  range  of  data  values . 

Figures  4a  and  4b  show  the  geographical  distribution  of  cloud  cover 
determined  from  observation  and  simulation,  respectively.  The  corresponding  gray 
shade  values  are  depicted  in  the  lower  right-hand  corner.  The  brighter  the  cloud 

picture,  the  larger  the  cloud  cover  is. 

The  major  distinctive  features  of  the  geographical  distribution  of  cloud 
cover  depicted  in  the  observation  and  simulation  include  the  following.  First, 
a  pronounced  bright  band  of  cloudiness  extends  along  the  ITCZ  from  150*  W  to  90* 

E.  This  cloud  band  splits  into  two  bands.  The  north  branch  spreads  into 
southern  Asia  and  the  northwestern  Pacific,  the  south  branch  expands  southeast 
to  the  southwestern  Pacific.  The  north  band  is  associated  with  the  seasonal 
monsoon,  while  the  south  band  is  related  to  the  Southern  Pacific  Convergence  Zone 
(SPCZ).  Although  the  ITCZ  cloud  band  is  less  pronounced  in  the  simulation,  the 
model  reasonably  predicts  the  bright  cloud  bands  extending  from  the  central 
Pacific  to  the  central  Atlantic,  and  into  central  Africa. 

Second,  areas  of  small  cloudiness,  i.e.,  a  dark  cloud  picture,  associated 
with  subtropical  highs  occur  both  north  and  south  of  the  ITCZ.  For  example,  the 
small  cloud  covers  over  southern  Africa,  Europe  and  Middle-East  Asia,  and  North 
and  South  America  are  clearly  identified  in  both  the  observation  and  the 
simulation.  The  oceanic  subtropical  highs  in  the  Atlantic  and  Pacific  Oceans  are 
also  evident  from  the  presence  of  small  cloudiness  over  the  regions. 

Third,  the  large  bright  cloud  bands  in  the  mid- latitudes  of  both 
hemispheres  are  associated  with  the  mid-latitude  storm  tracks.  In  both  the 
simulation  and  the  observation,  bright  cloud  bands  are  found  over  the  east  cost 
of  South  America,  the  southern  and  northern  Pacific,  and  the  Greenland  Sea.  The 
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cloudiness  over  the  east  cost  of  North  America  is  larger  in  the  simulation  than 
in  the  observation. 

Finally,  the  under-  and  overestimation  of  cloud  cover  is  shown  in  the 
3DNEPH  data  in  the  Arctic  and  Antarctic  regions,  respectively,  for  the  Northern 
Hemisphere  summer.  Comparing  Figs.  4a  and  4b,  it  is  evident  that  the  cloud 
picture  over  the  Antarctic  area  is  brighter  in  the  observation  than  in  the 
simulation,  and  vice  versa  for  the  Arctic  region. 

Figures  5a  and  5b  show  the  OLR  obtained  from  observation  and  simulation, 
respectively.  The  intensity  of  the  gray  shades,  given  in  the  lower  right-hand 
corner  of  each  figure,  is  defined  such  that  the  brighter  the  picture,  the  smaller 
the  radiative  flux  is.  There  are  two  distinct  features  on  these  maps.  First, 
the  well  defined  ITCZ  and  monsoon  circulation  are  characterized  by  bright 
regions,  i.e.,  relatively  small  values  of  OLR.  They  are  located  over  the 
equatorial  western  and  eastern  Pacific,  the  Atlantic  and  central  Africa,  as  well 
as  the  monsoon  area  over  India  and  Asia.  As  mentioned  previously,  the  cloudiness 
in  these  regions  is  underestimated  in  the  simulation.  As  would  be  expected,  the 
flux  values  in  these  regions  are  overestimated  by  the  model.  Second,  the  almost 
continuous  dark  belts  are  associated  with  the  subtropical  highs  located  to  the 
south  and  north  of  the  ITCZ.  In  these  regions,  large  OLR  values  are  due  to  the 
small  cloudiness  in  the  areas.  Small  flux  values  are  also  evident  in  the  mid¬ 
latitude  storm  tracks  in  the  Northern  Hemisphere.  However,  the  simulated  values 
of  OLR  in  the  Greenland  Sea,  near  70"  N  and  0®  E,  appear  to  be  too  small  in 
comparison  with  observation. 

Figures  6a  and  6b  illustrate  the  geographical  distribution  of  LWP  derived 
from  observation  and  simulation,  respectively.  The  observed  data  taken  from 
Prabhakara  and  Short  (1984)  are  displayed  in  terms  of  the  monthly  averaged  map 
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for  July,  1979.  In  Figs.  6a  and  6b,  a  large  amount  of  liquid  water  in  the  ITCZ 
is  evident.  However,  in  the  subtropical  highs  of  both  hemispheres,  liquid  water 
is  small.  In  the  southern  Pacific,  the  model  simulates  large  LWPs ,  which  are 
also  evident  in  the  observed  values.  In  both  the  simulation  and  the  observation, 
abundant  liquid  water  is  found  over  the  Indian  Ocean,  and  southeast  and  northeast 
Asia.  In  addition,  the  well-defined  observed  local  LWP  maxima  near  the  east 
coasts  of  North,  Central,  and  South  America  are  well  simulated  by  the  model. 
Local  LUP  minima  observed  west  of  Africa,  and  South  and  North  America  are  also 
shown  in  the  simulation. 

The  large  amount  of  liquid  water  associated  with  the  mid- latitude  storm 
tracks  near  50’  S(N)  -  60’  S(N)  is  not  shown  in  the  SMMR  observation  because  the 
sea  surfaces  are  largely  covered  by  ice.  In  addition,  the  model  simulates  large 
amounts  of  liquid  water  over  land  in  Central  America,  central  Africa,  India  and 
southeast  China.  These  results  are  not  reported  by  the  SMMR  observation  because 
of  the  inability  of  the  microwave  technique  to  function  over  land.  From  the 
above  comparisons,  it  appears  that  the  geographic  distribution  of  LWP  is  well 
simulated  in  the  present  model. 

Since  observed  cloud  IWC  has  not  been  available  on  the  global  scale,  we 
qualitatively  compare  the  model -predicted  cloud  IWC  with  values  derived  from 
aircraft  measured  data  in  cirrus,  as  well  as  with  other  model  results.  Based  on 
aircraft  measured  data  in  mid- latitude  cirrus  over  the  United  States,  Heymsfield 
and  Platt  (1984)  have  derived  the  ice  crystal  size  distribution  and  IWC  as 
functions  of  temperature.  Averaged  cloud  IWC  varies  from  0.001  g/m^  at 
temperatures  below  -40’C  to  about  0.02  g/m^  above  this  temperature.  Similar 
cloud  IWC  values  have  also  been  simulated  from  a  cirrus  cloud  model  by  Starr  and 
Cox  (.1985).  The  simulated  cloud  IWC,  averaged  over  July  3-4,  1979,  is  shown  in 
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Fig.  7.  The  magnitude  of  the  predicted  cloud  IWC  is  on  the  order  of  0.01  g/m^, 
which  is  consistent  with  the  aforementioned  model  simulation  and  aircraft 
observed  values  in  mid- latitude  cirrus  clouds. 

The  maxima  cloud  IWC  are  located  along  a  temperature  of  about  -15*C,  at 
which  the  maximum  rate  of  Bergeron-Findeison' s  process  occurs.  The  simulated 
cloud  IWC  decreases  from  about  0.01  g/m^  near  -15*C  to  0.001  g/n?  at  temperatures 
below  -40®C.  The  variations  of  cloud  IWC  with  respect  to  the  temperature 
obtained  from  the  model  are  similar  to  those  suggested  by  Heymsfield  and  Platt 
(1984).  It  is  noted  that  the  simulated  cloud  ice  vanishes  at  temperatures  warmer 
than  the  melting  point,  because  the  melting  process  converts  ice  crystals  into 
water  droplets.  The  simulated  cloud  IWC  appears  to  depend  more  on  temperature 
than  on  vertical  velocity,  which  is  most  significant  in  the  condensation  process. 
The  latitudinal  distribution  of  cloud  IWC  appears  to  be  realistic.  In  the 
wintertime  Antarctic  area,  all  clouds  contain  ice,  while  in  the  summertime  Arctic 
region,  only  middle  and  high  clouds  are  ice  clouds.  In  the  tropics,  high  clouds 
are  largely  composed  of  ice  crystals. 


5.  THE  ROLE  OF  RADIATIVE  HEATING  AND  ICE  PHASE  PROCESSES  IN  CLOUD  FORMATION 

In  order  to  examine  the  effects  of  radiative  heating  fields  and  ice  phase 
microphysics  on  the  large-scale  cloud  simulations,  a  set  of  numerical  experiments 
has  been  carried  out  with  and  without  the  inclusion  of  radiative  heating  and  ice 
phase  processes.  Below  is  a  brief  description  of  the  experimental  runs; 

•  CTRL:  Control  run  which  includes  both  radiation  and  ice -phase 

processes . 

•  EXPl:  As  in  CTRL,  but  without  radiative  heating. 

•  EXP2 :  As  in  CTRL,  but  without  the  ice-phase  processes. 


All  of  the  control  and  experimental  runs  have  the  same  initial  conditions  and 
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wind  fields.  The  results  are  averaged  over  the  same  period,  allowing 
intercomparisons  to  be  made  between  the  control  and  experimental  runs. 

a .  Control  R\m 

The  control  run  is  identical  to  the  experiment  described  in  Section  4b. 
In  that  section,  the  computational  results  were  illustrated  in  terms  of  the 
horizontal  cross  section.  In  this  section,  these  results  are  shown  in  terms  of 
the  meridional  cross  section  so  that  the  vertical  structure  of  cloud  fields  can 
be  easily  identified. 

Figures  8a  and  8b  show  the  zonally  averaged  cloud  cover  and  cloud  LUC 
predicted  from  the  model  CTRL.  Cloud  IWC  predicted  from  CTRL  has  been  shown  in 
Fig.  7.  Three  cloud  cover  maxima  are  seen  in  the  tropics  and  in  the  mid-high 
latitudes  of  both  hemispheres.  These  maxima  are  due  to  large  upward  motions  in 
these  regions.  Minima  cloud  covers  occur  in  the  regions  of  subtropical  highs, 
which  are  directly  related  to  downward  motions.  Because  of  the  strong  sinking 
motion  of  the  Southern  Hemisphere  Hadley  cell,  cloud  cover  in  the  subtropical 
high  of  the  Southern  Hemisphere  is  smaller  than  that  in  the  Northern  Hemisphere 
counterpart. 

The  pattern  of  simulated  cloud  LUC  in  Fig.  8b  is  consistent  with  that 
illustrated  in  Fig.  8a.  The  three  cloud  LUC  maxima,  located  in  the  tropics  and 
mid- latitude  storm  tracks,  correspond  to  large  cloudiness  in  these  regions.  The 
cloud  LUC  minima  associated  with  the  subtropical  highs  are  due  to  the  prevailing 
downward  motions  in  the  regions.  In  the  vertical,  large  cloud  LWCs  are 
associated  with  low  and  middle  clouds.  The  maximum  cloud  LWC  is  located  in  the 
tropics,  and  two  local  maxima  exist  in  mid- latitudes . 

b.  Experiment  on  Inclusion  of  Radiative  Heating 

Figures  11-13  show  the  results  from  the  experiments  with  and  without  the 
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inclusion  of  radiative  heating.  These  results  are  expressed  in  terms  of  the 
differences  between  two  experiments  (CTRL-EXPl)  . 

Solar  heating  rate  differences  are  illustrated  in  Fig.  9a.  Solar  radiation 
heats  the  atmosphere.  As  would  be  expected,  the  solar  heating  rates  in  July  of 
the  Northern  Hemisphere  are  much  larger  chan  those  of  the  Southern  Hemisphere 
counterpart.  There  is  no  solar  radiation  within  the  Antarctic  Circle  because  the 
sun  does  not  rise  above  the  horizon  in  this  region  during  July.  Two  solar 
heating  maxima  with  values  of  -2  K/day  are  located  at  the  tops  of  middle  clouds 
in  the  Northern  Hemisphere.  Large  solar  heating  rates  within  the  Arctic  Circle 
are  due  to  long  solar  days.  Large  solar  heating  rates  in  the  tropics  are  related 
to  small  solar  zenith  angles  during  the  daytime.  The  large  solar  heating  located 
in  the  upper  atmosphere  is  associated  with  strong  ozone  absorption. 

Figure  9b  shows  the  IR  heating  rate  differences.  The  most  distinct  feature 
in  this  figure  is  the  strong  IR  cooling  rates  that  occur  at  the  middle  cloud 
tops.  This  is  because  high  clouds  are  considered  to  be  nonblack  in  radiation 
calculations,  while  middle  clouds  are  assumed  to  be  blackbodies.  As  a  result, 
the  upward  IR  flux  emitted  by  middle  clouds  is  much  stronger  than  the  downward 
flux  emitted  by  high  clouds.  Consequently,  strong  IR  radiative  cooling  is 
generated  at  the  tops  of  middle  clouds. 

The  inclusion  of  IR  radiative  heating  basically  cools  the  atmosphere. 
However,  if  low  cloudiness  is  sufficiently  large  to  emit  strong  downward  IR 
fluxes  from  the  cloud  base,  weak  IR  heating  may  be  produced  beneath  the  cloud 
base.  This  explains  the  two  regions  of  weak  IR  heating  at  the  low  cloud  bottoms 
at  60*  S  and  in  the  tropics.  Two  other  weak  IR  heating  rates,  located  in  the 
upper  tropical  atmosphere  and  in  the  lower  atmosphere  in  the  Arctic,  are 
associated  with  the  areas  of  high  ozone  concentration  and  temperature  inversion 
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in  the  regions,  respectively. 

Figure  10a  shows  net  radiative  heating  rates.  The  inclusion  of  the 
radiative  processes  in  the  model  decreases  the  atmospheric  temperature.  As  would 
be  expected,  the  net  radiative  heating  shows  relatively  large  cooling  at  middle 
cloud  tops  and  small  cooling  at  cloud  bottoms.  However,  there  are  weak  heating 
rates  in  lower  tropical  atmospheres  and  in  the  Southern  Hemisphere  due  to  large 
low  cloud  cover  over  those  regions .  The  strongest  radiative  cooling  at  the 
middle  cloud  top,  located  within  the  Antarctic  Circle,  is  due  to  the  absence  of 
solar  radiation  in  the  winter  season.  Relatively  large  cooling  at  the  low  cloud 
level  within  the  Arctic  area  is  associated  with  the  temperature  inversion  in  the 
region,  while  the  large  cooling  near  the  surface  layers  in  the  equatorial  region 
is  due  to  the  contribution  of  water  vapor  continuum  absorption.  It  is  noted  that 
noticeable  net  radiative  heating  is  located  in  the  upper  atmosphere  in  the 
Northern  Hemisphere  because  of  the  ozone  solar  absorption  in  the  region.  The 
weak  heating  at  the  middle  cloud  level  within  the  Arctic  Circle  is  due  to  the 
absorption  of  solar  radiation  during  long  solar  days. 

Figure  10b  shows  the  temperature  pattern.  The  inclusion  of  radiative 
heating  decreases  the  atmospheric  temperature  which  decreases  more  significantly 
for  cloudy  conditions  than  clear  sky.  The  largest  temperature  reduction  occurs 
at  the  middle  cloud  top  near  66.5*  S,  where  the  sun  does  not  rise  much  above  the 
horizon  in  the  winter  season.  Temperature  decreases  are  relatively  small  in  the 
subtropical  regions  where  small  cloud  amounts  are  located.  These  temperature 
changes,  in  turn,  initiate  all  of  the  changes  in  the  cloud  fields  such  as  cloud 
cover,  cloud  LWC,  and  cloud  IWC. 

Figure  11a  shows  differences  in  cloud  cover  between  the  two  experiments. 
The  inclusion  of  radiative  heating  leads  to  an  increase  in  cloud  cover.  Large 
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increases  occur  at  the  middle  cloud  tops,  at  which  large  temperature  decreases 
are  located.  The  increases  in  low  cloudiness  are  relatively  small  due  to  small 
radiative  cooling  associated  with  low  clouds.  However,  there  are  two  large 
increases  of  low  cloud  cover.  One  occurs  in  the  tropics  where  a  large  amount  of 
water  vapor  exists,  while  the  other  is  in  the  Arctic  area  where  a  strong 
inversion  is  located. 

Figure  11b  illustrates  the  differences  in  total  cloud  water  content,  i.e., 
cloud  LWC  plus  IWC  between  CTRL  and  EXPl.  Cloud  water  increases  wherever  cloud 
cover  increases.  As  has  been  mentioned  previously,  the  decrease  of  cloud 
temperature  by  radiative  cooling  increases  the  condensational  rate  and  results 
in  the  increase  of  total  cloud  water.  The  large  increase  of  cloud  water  occurs 
at  the  tops  of  middle  and  low  clouds  in  the  tropics,  where  strong  radiative 
cooling  is  located. 

c.  Experiment  on  Ice  Phase  Processes 

Figures  12a  and  12b  show  the  results  from  the  experiments  with  and  without 
the  inclusion  of  ice -phase  processes.  The  differences  in  temperature  and  cloud 
cover  due  to  the  inclusion  of  ice-phase  processes  are  insignificant.  Figure  12a 
illustrates  the  differences  in  cloud  LWC  between  the  cases  with  and  without  the 
inclusion  of  ice -phase  processes.  This  figure  shows  little  change  in  cloud  LWC 
in  the  tropical  low  and  middle  clouds  that  contain  primarily  water  droplets. 
Marked  reductions  of  cloud  LWC  with  ice  processes  included  occur  in  low,  middle, 
and  high  clouds  within  the  Antarctic  Circle,  high  clouds  in  the  tropics,  and 
middle  and  high  clouds  within  the  Arctic  Circle.  In  these  regions  clouds  are 
largely  composed  of  ice  crystals,  which  deplete  the  cloud  LWC  by  virtue  of 
Bergeron-Findeison' s  process.  Since  there  is  no  cloud  IWC  in  EXP2,  the 
differences  in  cloud  IWC  are  identical  to  cloud  IWC  shown  in  CTRL  (see  Fig.  7  ). 
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Figure  12b  shows  the  differences  in  total  cloud  water  content.  The 
increases  of  total  cloud  water  when  the  ice  phase  is  included  indicate  that  the 
increase  of  ice  crystals  in  clouds  is  larger  than  the  reduction  of  cloud  LWC. 
It  follows  that  ice  clouds  in  these  regions  are  not  entirely  controlled  by 
Bergeron-Findeison' s  process. 

In  order  to  further  examine  the  relationship  between  autoconversion  and 
Bergeron-Findeison’ s  process  on  cloud  formation,  we  have  undertaken  two  more 
experimental  iruns  with  the  autoconversion  rate  of  Cq  -  10"^  sec'\  which  is  five 
times  smaller  than  what  was  used  in  the  CTRL.  A  smaller  autoconversion  rate 
leads  to  the  production  of  more  cloud  LWC,  which  in  turn  generates  more  cloud  IWC 
through  Bergeron-Findeison' s  process.  In  the  following  two  experimental  runs, 
the  first  run  is  the  s-ime  as  CTRL  with  the  smaller  value  of  the  autoconversion 
rate,  i.e.,  Cq  -  ,  the  second  run  is  the  same  as  the  first  run  but  without 
the  ice  phase  processes. 

Figure  13a  shows  a  pattern  very  similar  to  that  in  Fig.  12a.  Due  to  the 
smal'^er  value  of  the  autoconversion  rate,  the  differences  in  cloud  LWC  are  about 
five  times  larger  in  Fig.  13a  than  they  are  in  Fig.  12a.  The  large  differences 
in  cloud  LWC  indicate  that  a  strong  Bergeron-Findeison  process  takes  place  when 
there  is  more  cloud  LWC.  However,  the  values  of  cloud  IWC  in  Fig.  13b  are  close 
to  those  in  Fig.  7,  even  though  Bergeron-Findeison' s  process  is  much  stronger  in 
the  former  than  in  the  latter.  This  indicates  that  the  downward  ice  flux  is 
larger  in  the  former  than  in  the  latter. 

Figure  13c  shows  the  differences  in  total  cloud  water  content  for  the 
smaller  Cg.  The  total  cloud  water  content  decreases  in  the  region  where  ice 
clouds  are  present.  The  reduction  of  the  total  cloud  water  content  is  mainly 
controlled  either  by  autoconversion  that  converts  cloud  LWC  into  precipitation 
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or  by  gravitational  settling  that  removes  ice  crystals  from  the  cloud  layer.  The 
decrease  of  total  cloud  water  content  when  the  ice  phase  processes  are  included 
indicates  that  gravitational  settling  is  an  efficient  mechanism  for  reducing  the 
cloud  IWC. 

6.  CONCLUSIONS 

A  time -dependent,  three-dimensional,  large-scale  cloud-moisture  model  has 
been  developed  for  the  prediction  of  temperature,  cloud  cover,  c?  LWC ,  cloud 
IWC,  and  precipitation.  Partial  cloudiness  is  allowed  to  form  when  the  large- 
scale  relative  humidity  is  less  than  100%.  The  cloud  microphysical  processes 
simulated  in  the  model  include  both  liquic.  and  ice  phases,  which  are  physically 
included  for  the  first  time  in  a  large-scale  cloud  model  to  simulate  cirrus 
clouds . 

A  96-hour  model  simulation,  with  the  initial  conditions  taken  at  12Z  on 
July  1,  1979  has  been  carried  out  to  assess  the  performance  of  the  large-scale 
cloud  model.  The  comparisons  between  the  model  simulation  and  the  satellite 
observations  indicate  the  large-scale  cloud  model  is  capable  of  realistically 
simulating  zonal  means  and  geographical  distributions  of  cloud  fields,  including 
cloud  cover,  OLR  and  net  solar  flux  at  TOA,  and  cloud  LWC.  Differences  in  the 
zonal  mean  total  cloud  cover  between  the  simulation  and  the  3DNEPH  data  are 
generally  within  10%.  Generally,  the  simulated  zonal  mean  OLRs  are  within  10 
W/m^  of  the  ERB  results.  The  zonal  mean  values  of  simulated  cloud  LWP  compare 
quite  well  with  the  SMMR  observations,  with  differences  being  less  than  50  g/m^, 
which  is  the  uncertainty  of  the  SMMR  data.  However,  it  is  clearly  indicated  that 
the  model  underestimates  cloud  cover  and  cloud  LWC  and  overestimates  the  OLR  in 
the  tropics,  where  clouds  are  predominantly  convective  clouds,  which  are  not 
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simulated  in  this  model.  In  the  geographical  distributions  of  cloud  fields,  the 
large-scale  cloud  model  has  simulated  the  large  cloud  bands  in  the  ITCZ,  monsoon 
areas,  Southern  Pacific  Convergence  Zone,  and  the  mid-latitude  storm  tracks, 
although  the  ITCZ  cloud  band  is  less  pronounced  in  the  simulation.  At  the 
subtropical  highs,  areas  of  small  cloudiness  over  southern  Africa  and  northern 
Africa,  and  North  and  South  America  are  clearly  identified  in  both  of  the 
observations  and  the  model  simulation.  The  magnitude  of  the  cloud  IWC  in  the 
simulation  is  about  0.01  g/m^,  which  is  consistent  with  the  aircraft  measurements 
in  cirrus  clouds,  as  well  as  with  other  model  results.  The  simulated  cloud  IWC 
appears  to  depend  more  on  temperature  than  on  vertical  velocity,  which  is 
critical  to  the  condensational  process. 

The  model-simulated  cloud  fields  are  strongly  linked  to  large-scale 
circulations.  Inspecting  the  model-predicted  cloud  fields,  such  as  cloud  cover, 
OLR,  and  cloud  LWC,  we  find  that  these  cloud  properties  are  closely  related  to 
one  another.  Large  amounts  of  cloud  cover  coincide  with  abundant  cloud  LWC, 
which  generates  large  precipitation  rates  by  virtue  of  the  autoconversion  process 
simulated  in  the  model.  In  summary,  the  atmospheric  moisture  and  cloud  fields 
are  linked  together  by  the  cloud  microphysical  processes  simulated  in  the  cloud 
model.  More  importantly,  the  moisture,  cloud  fields  and  earth's  radiation  budget 
are  physically  related  to  large-scale  thermal  and  dynamical  structures,  such  as 
temperature  and  vertical  velocity. 

Sensitivity  studies  have  been  performed  to  examine  the  effects  of  radiative 
heating  fields  and  ice-phase  cloud  microphysics  on  large-scale  cloud  formation. 
The  inclusion  of  radiative  processes  in  the  cloud  formation  model  significantly 
decreases  the  cloud  temperature,  with  the  strongest  cooling  occurring  at  the 
middle  cloud  top.  However,  there  may  be  weak  heating  located  in  a  low  cloud  if 
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cloud  cover  is  sufficiently  large  to  produce  a  strong  downward  flux.  The 
decrease  of  temperature  by  radiative  cooling  increases  the  atmospheric  relative 
humidity  and  condensation  in  clouds.  Af  a  result,  cloud  cover  and  total  cloud 
water,  including  cloud  LWC  and  IWC,  increase  due  to  the  inclusion  of  radiative 
processes . 

The  inclusion  of  ice-phase  processes  has  a  greater  effect  on  the  cloud  LWC 
and  IWC  than  on  the  temperature  and  cloud  cover  fields.  If  the  cloud  LWC  is 
large  enough  to  undertake  an  effective  Bergeron-Findeison' s  process,  the  total 
cloud  water  content  decreases,  with  large  decreases  located  at  the  regions  where 
the  presence  of  ice  clouds  is  most  pronounced.  The  reduction  of  the  total  cloud 
water  content  is  mainly  controlled  either  by  autoconversion,  which  converts  cloud 
LWC  into  precipitation,  or  by  gravitational  settling,  which  removes  ice  crystals 
from  the  cloud  layer.  The  decrease  of  total  cloud  water  content  when  the  ice- 
phase  processes  are  included  indicates  that  gravitational  settling  is  an 
efficient  mechanism  for  reducing  the  cloud  IWC. 
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FIGURE  CAPTIONS 


Figure  1.  Empirical  relations  between  the  low-,  middle-,  and  high- cloud  amount 
and  relative  humidity  (Smagorinsky,  1960).  The  values  of  hg  for 
low,  middle,  and  high  clouds  are  0.6,  0.35,  and  0.26,  respectively. 

Figure  2.  Schematic  display  of  the  cloud  microphysical  processes  simulated  in 
the  model. 

Figure  3.  (a)  Comparison  of  zonal-mean  total  cloudiness  between  3DNEPH  and 

model  simulated  results  averaged  over  the  time  period,  July  3-4; 
(b)  Comparison  of  the  zonally  averaged  OLR  computed  from  the  model 
and  derived  from  ERfi  data.  The  model  results  and  ERB  data  are 
averaged  over  the  time  period,  July  3-4;  (c)  Comparison  of  the 
zonal  mean  LWP  in  the  simulation  and  from  SMMR  data.  The  SMMR 
results  are  taken  from  Prabhakara  and  Short  (1984)  for  the  monthly 
average  of  June,  1979.  Only  the  LWPs  that  are  zonally  averaged  over 
the  oceans  are  shown  in  the  figure  for  both  simulation  and  SMMR 
data. 

Figure  4.  The  geographical  distribution  of  total  cloud  cover  averaged  over  the 
period,  July  3-4  (a)  derived  from  3DNEPH  and  (b)  simulated  by  the 
model.  The  scale  on  the  right  of  each  figure  shows  the  shading  for 
the  different  cloud  amount  intervals.  The  brighter  the  cloud 
picture,  the  larger  the  cloud  cover  is. 

Figure  5.  The  geographical  distribution  of  OLR  averaged  over  the  time  period, 
July  3-4  in  V/m  (a)  derived  from  ERB  data  and  (b)  simulated  by  the 
model.  The  scale  on  the  right  of  each  figure  shows  the  shading  for 
the  different  flux  intervals. 

Figure  6.  The  geographical  distribution  of  LWP  in  10  g/m^  (a)  derived  from 
-  SMMR  data  for  the  monthly  average  of  July,  1979  (Prabhakara  and 
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Short,  1984)  and  (b)  simulated  by  the  model  for  two  days  average 
covering  July  3-4. 

The  model  simulations  for  zonal -mean  cloud  IWC  with  a  contour 
interval  of  2.5  x  10“^  g/®^- 

The  model  predicted  zonal  mean  fields  of  (a)  cloud  cover,  and  (b) 
cloud  LWC.  The  model  results  are  averaged  over  the  time  period, 
July  3-4. 

Zonally  averaged  differences  in  solar  and  IR  radiative  heating 
rates.  The  differences  are  obtained  by  subtracting  the  results  in 
EXPl  from  those  in  CTRL.  The  contour  intervals  for  solar  and  IR 
heating  rates  are  0.5  and  2  K/day,  respectively. 

As  in  Fig.  9,  but  for  differences  in  net  heating  rate  and 
temperature . 

As  in  Fig.  9,  but  for  differences  in  cloud  cover  and  total  cloud 
water  with  contour  intervals  of  10%  and  2  x  10"^  g/®^,  respectively. 
Zonally  averaged  differences  in  cloud  LWC  and  total  cloud  water. 
The  differences  are  obtained  by  subtracting  the  results  in  EXP2  from 
those  in  CTRL.  The  contour  intervals  for  cloud  LWC  and  total  cloud 
water  are  2  x  10"^  g/m  and  10'^  g/®^,  respectively. 

Zonally  averaged  differences  in  (a)  cloud  LWC,  (b)  cloud  IWC,  and 
(c)  total  cloud  water.  The  differences  are  obtained  by  subtracting 
the  results  without  ice  phase  from  those  with  ice  phase  in  the  case 
of  a  small  autoconversion  rate.  The  contour  intervals  for  cloud  LWC 
and  total  cloud  water  are  2  x  10'^  g/m^  and  10'^  g/®^,  respectively. 
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Summary 

Parameterization  programs  for  cloud  microphysics  and  ra¬ 
diative  transfer  involving  ice  clouds  have  been  developed  in 
tenm.s  of  the  mean  effective  size  and  ice  vsater  path  The  mean 
effective  size  appears  to  be  adequate  m  representing  the  ice 
crystal  size  distribution  for  radiative  parameterizations.  For 
a  given  ice  water  path,  smaller  mean  effective  sizes  reflect 
more  solar  radiation,  emit  more  IR  radiation  and  enhance 
net  radiative  heating  cooling  at  the  cloud  top  and  bottom 
than  larger  sizes.  The  presence  of  small  ice  crystals  may  gen¬ 
erate  steeper  lapse  rates  in  clouds  A  3-D  global  cloud  model 
that  prescribes  the  horizontal  wind  fields  in  a  24  hour  period 
is  used  to  investigate  the  sensitivity  of  the  mean  effective  size 
of  ice  crystals  on  the  simulation  of  radiative  heating,  tem¬ 
perature.  cloud  cover  and  ice  water  content  A  variation  in 
the  mean  effective  size  from  75  to  50  pm  in  a  24  hour  pre¬ 
diction  simulation  generates  more  cooling  above  the  high 
cloud  top  and  a  decrease  of  temperature  These  results  lead 
to  an  increase  of  high  cloud  cover  in  some  latitudes  by  as 
much  as  4°u  and,  at  the  same  time,  a  dec. case  of  middle 
cloud  cover  bv  3  4%  in  latitudes  between  60”S  and  60'N. 

1.  Introduction 

In  recognition  of  the  importance  of  the  role  of 
clouds  in  weather  and  climate  processes,  numerous 
efforts  have  been  made  to  include  a  more  realistic 
and  physically  based  parameterization  for  cloud 
parameters  in  general  circulation  models  (Sundqv- 
ist.  1988;  Liou  and  Zheng.  1984;  Heymsfield  and 
Donner.  1990;  Smith,  1990),  To  understand  the 
effects  and  feedbacks  of  clouds  on  dynamic  and 


thermodynamic  systems,  the  microphysical  pro¬ 
cesses  that  govern  the  formation  and  dissipation 
of  clouds  must  be  adequately  incorporated  in  at¬ 
mospheric  models.  Furthermore,  the  radiative 
properties  of  clouds  are  determined  by  the  cloud 
thermodynamic  phase  (ice  or  water),  particle  size 
distribution  and  geometric  configuration.  It  ap¬ 
pears  to  be  impractical  to  include  individual  cloud 
particles  in  cloud  models  and  radiation  calcula¬ 
tions.  We  must  rely  on  some  aspect  of  the  cloud 
particle  size  distribution  in  the  performance  of 
radiation  parameterizations  in  numerical  models. 

In  this  paper,  we  e.vplore  the  means  by  which 
cloud  microphysics  and  radiative  transfer  can  be 
effectively  and  physically  parameterized  for  in¬ 
corporation  in  large-scale  cloud  and  atmospheric 
models.  We  shall  confine  our  presentation  to  high 
clouds  that  contain  ice  particles.  In  Sections  2  and 
3,  we  present  parameterization  programs  for  cloud 
microphysics  and  radiative  transfer  in  clouds,  re¬ 
spectively .  The  mean  effective  size  is  proposed  to 
represent  ice  crystal  size  distribution  in  radiation 
parameterizations.  In  Section  4,  we  investigate  the 
sensitivity  of  mean  effective  size  to  the  simulation 
of  radiative  heating,  temperature,  cloud  cover  and 
ice  water  content  (IWC)  using  the  global  cloud 
model  developed  b\  Lee  et  al.  (1990).  Finally, 
conclusions  are  given  in  Section  5. 
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2.  Parameterization  of  Cloud  Microph\sics 

It  is  not  computationally  feasible  to  account  for 
individual  cloud  particles  in  radiative  transfer  cal¬ 
culations.  Nor  is  it  practical  to  incorporate  a  pro¬ 
gram  that  takes  into  account  the  growth  of  in¬ 
dividual  water  droplets  and  ice  crystals  in  cloud 
models.  However,  some  measure  of  the  particle 
size  distribution  may  be  used.  From  the  perspec¬ 
tive  of  radiation  calculations,  particles  scatter  an 
amount  of  light  proportionate  to  their  cross-sec¬ 
tion  area.  The  cross-section  of  a  nonspherical  par¬ 
ticle  is  proportional  to  the  ma.'timum  dimension 
(or  the  length.  L)  multiplied  by  the  width,  D.  Thus, 
we  may  define  a  mean  effective  width  (size)  for 
he.xagonal  ice  crystals  in  the  form 

D,  =  J  D  LDn(L)dLi  J  LDn{L)dL  (I) 

^Tnin  ^min 

where  and  denote  the  limits  of  the  ice 
crystal  length.  From  laboratory  and  aircraft  ob¬ 
servations  (Auer  and  Veal,  1970;  Heymsfield, 
1972;  Hobb  et  al..  1974),  the  width  of  an  ice  par¬ 
ticle  can  be  related  to  the  length  via  the  following 
equation;  D  =  ayV'\  where  a\  and  6,  are  certain 
coefficients  that  are  dependent  on  the  particle 
shape. 

The  ice  water  content  is  defined  by 
I\VC=  I  m(L)n(L)dL.  (2) 


where  the  mass  of  an  individual  ice  crystal  may 
be  related  to  the  length  (Heymsfield.  1972);  ni  = 
with  u;  and  being  certain  coefficients. 
From  Eqs.  ( 1 )  and  (2),  IWC  and  D,  may  be  related 
through  the  ice  crystal  size  distribution.  The  ver¬ 
tical  IWC  or  ice  water  path  (IWP)  for  a  given 
thickness,  Jr,  is 

IWP  =  IWC  Jr.  (3) 

The  other  cloud  physics  parameter  that  is  critically 
dependent  on  ice  crystal  size  distribution  is  the 
bulk  terminal  velocity,  which  is  defined  by 

IV  =  n(L)  m(L)v{L)dL!  n{L)m(L)dL, 

Lrtwn  £.mia 

(4) 


where  i- iv  the  terminal  veloeitv  for  an  indi\idual 
ice  crystal,  which  can  he  determined  by  the  crvstal 
length  based  on  laboratory  results  (Javaweera  and 
Cottis,  1969;  Heymsfield  1972).  Thus  we  have:  v 
=  u,  L!".  with  u,  and  6-,  being  certain  coefficients 
(Starr  and  Cox,  1985). 

The  ice  crystal  size  distribution  for  midlatitude 
cirrus  clouds  has  been  measured  by  Heymsfield 
(1977)  using  an  optical  probe.  Based  on  the  mea¬ 
sured  data,  Heymsfield  and  Platt  (1984)  grouped 
the  size  distributions  according  to  temperatures 
ranging  from  -  20  to  -  60°C  and  proposed  a  par¬ 
ameterized  equation  in  the  form 

n{L)  =  AL^.  (5) 

where  both  coefficients  A  and  B  are  dependent  on 
temperature.  The  minimum  length  that  an  optical 
probe  can  measure  is  about  20  (am.  It  is  probable 
that  ice  crystals  smaller  than  20  pm  could  be 
missed  by  this  measurement  technique. 

In  view  of  the  preceding  discussion,  the  ice  crys¬ 
tal  width,  mass,  terminal  velocity  and  size  distri¬ 
bution  may  be  parameterized  in  terms  of  the 
length.  Using  Z)„,„  =  20  pm,  the  bulk  terminal 
velocity  and  mean  effective  width  can  be  expressed 
as  a  function  of  the  IWC,  as  shown  in  Fig.  1.  In 
the  calculations,  we  use  a  temperature  range  of 
( -  30 —  35'C)  to  obtain  the  coefficients  A  and  B. 
Both  Df  and  increase  with  increasing  IWC.  To 
investigate  the  effect  of  small  ice  crystals  on  these 
relationships,  in  Eq.  (1)  is  set  at  10  pm  and 
it  is  assumed  that  Eq.  (4)  is  valid  for  the  ice  crystal 
size  distribution  that  includes  small  ice  crystals 
with  sizes  from  10  to  20  pm.  The  bulk  terminal 
velocity  ranges  from  about  1  to  10  m  s  for  the 
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Fig,  I  Terminal  velocity  and  mean  elTecIive  size  as  functions 
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IWC'  ftoiii  0  1  to  I  c  Hi'.  For  IW  C's  snrallci  ihitii 
(i-lg  in',  the  bulk  terminal  veloeii)  is  less  than 
1  m  s.  The  addition  of  ice  crystals  smaller  than 
20|ani  does  not  affect  the  bulk  terminal  velocity 
significantb .  Based  on  our  present  parameteri- 
zation.s  for  cloud  microphysics,  the  mean  effecltve 
size  can  also  be  theoretically  computed  in  terms 
of  IW’C,  For  IWCs  larger  than  0.1  g  m',  the  ad¬ 
dition  of  10  20pm  sized  ice  crystals  results  in  no 
noticeable  changes  in  mean  effectne  sizes.  For  an 
IW'C  of  0.01  g  m’.  the  computed  mean  effective 
size  IS  about  50  pm.  With  additional  small  ice  crys¬ 
tals.  this  size  reduces*  to  about  40  pm.  It  should  be 
emphasized  that  the  relationship  between  the 
mean  effective  size  and  IWC  that  has  been  de¬ 
veloped  is  purely  theoretical  and  requires  verifi¬ 
cation  from  laboratory  and  aircraft  experiments. 

Light  scattering  and  absorption  calculations  re¬ 
quire  the  ice  crystal  size  distribution.  However,  for 
radiative  heating  and  llux  calculations,  additional 
information  concerning  the  cloud  thickness  and 
I\H'C  is  needed.  It  is  conventional  to  use  the  optical 
depth  to  represent  the  optical  properties  of  a  me¬ 
dia.  The  optical  depth,  r.  is  the  product  of  the 
extinction  coefficient  and  thickness.  The  extinc¬ 
tion  coefficient  is  the  product  of  the  average  ex¬ 
tinction  cross  section,  o,..  and  the  number  density 
of  cloud  particles.  In  the  limits  of  geometric  optics, 
and  assuming  that  ice  crystals  are  randomly  ori¬ 
ented  in  space,  we  have  <7,  =  4Z)  +  1)  2 

(Takano  and  Liou.  1989a).  This  expression  is  de¬ 
rived  based  on  the  optical  theorem  that  the  ex¬ 
tinction  cross  section  of  a  large  particle  is  twice 
its  geometric  cross  section.  In  this  case  we  find 

r  5  (o -t- h  D,)  IWP.  (6) 

w  here  a  and  b  are  certain  constants.  For  radiative 
transfer  calculations  associated  with  a  cloud 
model,  the  two  parameters  defined  in  Eqs.  ( 1 )  and 
(2)  appear  to  be  appropriate  and  sufficient. 

3.  Parameterization  of  Radiative  Transfer  in  Clouds 

For  the  computatlon^  of  fluxes  and  heating  rates 
within  cloud  layers,  wc  must  perform  spectral  in¬ 
tegration  covering  the  entire  solar  and  IR  spectra. 
These  spectra  may  be  divided  into  a  number  of 
spectral  intervals  according  to  the  location  of  the 
gaseous  absorption  band''  involving  cloud  parti¬ 
cles  and  gases,  primarily  water  vapor.  In  order  to 
resolve  the  variation  in  the  refractive  index  of  ice 
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and  to  account  for  the  gascoio  absorption,  six  and 
12  bands  are  selected  for  solar  and  thermal  IR 
regions,  respectivelv .  The  radiative  transfer  meth¬ 
odology  used  is  the  delta-four-stream  approxi¬ 
mation  developed  by  Liou  et  al.  (1988).  This  ap¬ 
proximation  can  achieve  relative  accuracy  within 
about  5‘'o  for  all  atmospheric  conditions.  The  so¬ 
lution  of  this  approximation  is  in  analytic  form 
so  that  the  computer  time  involved  is  minimal, 
and  the  method  can  be  readily  applied  to  inhom¬ 
ogeneous  media  by  matching  the  continuity  re¬ 
quirement  of  diffuse  intensities.  To  apply  this  ap¬ 
proximation.  four  expansion  coefficients  of  the 
phase  function  are  required. 

The  scattering  and  absorption  properties  of 
hexagonal  ice  crystals  whose  size  parameters  are 
greater  than  .30  w  ere  computed  from  the  geometric 
ray-tracing  technique  developed  by  Takano  and 
Liou  (1989a.  b).  For  size  parameters  that  are  less 
than  30.  wc  used  the  Mie-type  solution  for  sphe¬ 
roids  dev  eloped  by  Asano  and  Sato  ( 1 980).  Single¬ 
scattering  computations  were  made  for  a  number 
of  ice  crystal  size  distributions  reported  by  Heyms- 
field  and  Platt  (1984)  and  recently  derived  from 
the  FIRE  cirrus  experiments  (Heymsfield.  private 
communication). 

As  pointed  out  previously,  the  optical  depth 
may  be  expressed  in  terms  of  D,  and  IWP.  in  the 
form  given  in  Eq.  (6).  in  the  limits  of  geometric 
optics  for  scattering  processes.  The  following  gen¬ 
eral  parameterization  form  has  been  derived  for 
the  spectral  optica!  depth: 

T,  =  (^  ,  A'')  •  IWP.  (7) 

w  here  a,,  is  the  know  n  coefficient,  and  the  subscript 

/  (  =  1.2 . 18)  is  the  index  for  the  spectral  band. 

For  solar  wavelengths,  the  extinction  cross  sec¬ 
tions  would  be  the  same  by  virtue  of  the  geometric 
optical  theorem.  Thus,  for  a  given  thickness,  zJr. 
the  optical  depth  is  a  constant  in  the  solar  region 
based  on  the  geometric  ray  tracing  method  for 
hexagonal  ice  crystals. 

The  single-scattering  albedo  is  defined  as:  w  = 

I  -  (7^  (7, .  where  and  represent  the  absorption 
and  extinction  cross  sections,  respectively.  In  the 
lirnits  of  geometric  optics,  we  have  a,  =  3  (D  2)* 
(■y^3  4L  D)  2.  as  pointed  out  previously  The 
absorption  cross  section  is  generally  dependent  on 
particle  volume,  which  is:  1=  3-y'3  D'  Z.  8  in  the 
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case  of  a  hexagonal  particle.  Since  L  is  related  to 
D  based  on  laboratory  and  aircraft  observations, 
a)  should  be  a  function  of  D.  Thus,  it  is  appro¬ 
priate  to  represent  to  in  terms  of  a  poKnomial 
function  in  the  form 


} 


i^) 


where  h„  is  the  known  coefficient  through  nu¬ 
merical  fittings.  The  third  order  polynomial  is  ac¬ 
curate  within  about  3%. 

The  phase  function  should  also  be  a  function 
of  ice  crystal  size  distribution,  viz.,  the  mean  ef¬ 
fective  size.  It  follows  that  the  expansion  coeffi¬ 
cients  of  the  phase  function  may  be  written 


(■O  I  j  X  ^  n  l  j  ' 

«  =  () 


(9) 


where  the  index  I  =  1,2, 3,4  for  the  delta-four- 
stream  approximation,  oii  =  3g,  and  g  is  the 
asymmetry  factor.  Again,  we  find  that  the  third 
order  polynomial  expansion  is  sufficient  to  pro¬ 
vide  accuracy  within  3®o.  For  solar  wavelengths, 
the  phase  functions  for  hexagonal  ice  crystals  show 
pronounced  22  and  46°  halo  maxima.  For  this 
reason,  a  simple  representation  using  the  asym¬ 
metry  factor  for  the  transfer  of  solar  radiation  in 
ice  clouds  in  inadequate.  However,  in  the  thermal 
IR  wavelengths,  halo  features  are  largely  sup¬ 
pressed  due  to  absorption.  Thus,  we  may  use  the 
asymmetry  factor  to  represent  the  phase  function 
through  the  Henyey-Greenstein  function: 

m,  =  (2/+l)/  (10) 

Finally,  to  incorporate  the  forward  peak  contri¬ 
bution  in  multiple  scattering,  we  may  use  the  sim¬ 
ilarity  principle  for  radiative  transfer  to  adjust  the 
optical  depth,  single-scattering  albedo  and  phase 
function  in  the  forms  (Liou  et  al.,  1988) 


r'  =  r(l  -  fib).  (11a) 

lb'  =  {\  -  -  fib).  (11b) 

[fij,-/(2/+  D]  (1  -  A  (11c) 

where  the  fraction  of  scattered  energy  residing  in 
the  forward  peak,/  =  f'jj  4.  The  wavelength  index 
/  has  been  omitted  in  these  equations. 

Using  the  preceding  parameterizations  for 
cloud  physics  and  radiative  transfer.  Fig.  2  shows 
solar  reflectance  and  abNorptance  and  IR  emitt- 


ance  lor  cmissi\it\)  as  functions  of  D.  and  IWP. 
.\  cosine  of  the  '.olar  zenith  angle  of  0.5  is  used 
in  solar  radiative  transfer  calculations.  For  a  given 
D...  solar  reflectance  increases  with  increasing 
IWP.  But  for  a  given  IWP,  a  cloud  with  a  smaller 
Z),.  reflects  more  solar  radiation  because  of  the 
larger  effective  cross  section  area.  For  example, 
for  an  IWP  of  100 g  m  solar  reflectances  of 


g. 


F  ig  2  Sol.ir  rollccijncc  .ind  .ibsorptancc  and  IR  cnattancc 
(f.r  crni^snu\  )  as  fanenons  of  the  mean  elTective  M/e  .ind  lec 
A  a  ter  path 
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45 


Heating  rate  (k/hr)  Heating  rate  (k/hr)  Heating  rate  (k/hr)  Heating  rate  (k/hr) 


Fig  .V  Net  radiaiiNC  heating  rales  as  funeiions  of  heiglit  and  cosines  of  the  solar  zenith  angle.  //„.  for  mean  effective  sizes  of 
50  and  75  nm 


about  70,  55,  44  and  30“o  are  shown  for  D^.  of  25, 
50,  75  and  100  |im,  respectively.  For  solar  ab- 
sorptance,  small  ice  crystals  also  absorb  more  solar 
radiation  for  IWPs  up  to  75  g  m  ~  r  However,  ab¬ 
sorption  of  solar  radiation  also  depends  on  for¬ 
ward  scattering  of  cloud  particles.  Larger  particles 
would  have  stronger  forward  scattering  and  ab¬ 
sorb  more  solar  radiation  when  the  lU'P  is  larger 
than  about  75g  m“c  The  IR  emissivity  is  signif¬ 
icantly  dependent  on  both  D,,  and  IWP.  For  a 
given  IWP  of  20g  m“%  emissivity  differs  by  as 
much  as  40%  for  of  25  and  100  pm. 

The  effects  of  the  ice  crystal  mean  effective  size 
on  net  radiative  heating  are  illustrated  in  Fig.  3. 
In  the  calculations,  we  use  an  IWP  of  26g  m~*, 
which  corresponds  to  an  IWC  of  0.013g  m“  ’  for 
a  typical  cirrostratus  with  a  thickness  of  2  km.  This 
cloud  is  inserted  in  the  standard  atmospheric  tem¬ 
perature  and  humidity  profiles.  A  solar  constant 
of  1365  W  m~^  and  a  surface  albedo  of  0.1  are 
used  in  the  solar  radiation  calculation.  The  optical 
depths  for  D,.  of  50  and  75  pm  are  approximately 
2  and  1,  respectively.  Net  radiative  heating  rates 
that  correspond  to  four  cosines  of  the  solar  zenith 
angle  of  0,  0.167,  0.5  and  0.833  are  displayed  in 
Fig.  3.  U  hen  /z„  =  0,  only  IR  heating  cooling  takes 
place.  In  this  case,  a  smaller  of  50  pm  produces 
heating  cooling  rates  of  about  0.6-0.7'C  h  '  '  close 
to  the  bottom  and  top  of  the  cloud  layer.  For  a 
D,  of  75  pm.  the  heating  cooling  rates  at  the  cloud 
bottom  and  top  reduce  by  as  much  as  a  factor  of 
two.  This  reduction  is  due  to  the  fact  that  a  cloud 


having  smaller  Z),,s  emits  more  IR  radiation  at  the 
colder  cloud  top.  The  cloud  absorbs,  at  the  same 
time,  more  IR  radiation  emitted  from  the  warmer 
surface  and  atmosphere  below,  resulting  in  strong 
heating  at  the  cloud  base  region.  The  solar  heating 
rate  increases  with  decreasing  D^.  and  occurs  pri¬ 
marily  at  the  cloud  top.  where  IR  cooling  is  much 
more  pronounced.  For  =  0.833,  the  solar  heat¬ 
ing  rate  is  a  factor  of  about  three  smaller  than  the 
IR  cooling  rate  at  the  cloud  top.  It  is  evident  from 
this  illustration  that  the  cloud  field  is  largely  con¬ 
trolled  by  IR  heating/cooling  and  that  the  presence 
of  small  ice  crystals  may  generate  a  steeper  lapse 
rate  within  clouds. 

4.  The  Role  of  the  Mean  Effective  Cloud  Particle 
Size  in  Large-Scale  Cloud  Processes 

To  investigate  the  sensitivity  of  cloud  microphys¬ 
ics  and  radiative  transfer  paranieterizations  to 
large-scale  cloud  processes,  the  3-D  global  cloud 
model  developed  by  Lee  et  al.  (1990)  is  used.  This 
cloud  model  consists  of  thermodynamic  equations 
for  the  prediction  of  water  vapor,  cloud  liquid 
water  content  (LWC)  and  IWC,  precipitation  and 
temperature.  However,  the  wind  field  is  prescribed 
using  the  results  from  a  general  circulation  model 
(GCM).  The  model  has  a  12-layer  stretch  vertical 
coordinate,  w  hich  has  a  fine  resolution  in  the  lower 
levels  Cloud  covers  are  computed  from  the  thresh¬ 
old  method  To  compare  with  the  observed  cloud 
fields  and  to  perform  radiative  transfer  calcula- 
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Huns,  ihc  model  eeneraied  cloud  parameters  are 
strapped  into  lou,  middle  and  high  clouds.  Tiie 
cloud  micropln steal  processes  considered  in  the 
model  are  evaporation,  condensation,  autocon¬ 
version,  sublimation  deposition,  homogeneous 
nucleation.  Bergeron  processes,  and  gravitational 
setting  for  ice  crvstals.  .^s  described  in  Section  2. 
the  gravitational  setting  is  parameterized  in  terms 
of  the  bulk  terminal  velocity,  which  is  critical  to 
the  simulation  of  IVVC.  Radiative  transfer  para- 
meterizations  for  clouds  follow  those  presented  in 
Section  ?).  In  particular,  two  cloud  parameters,  i.e.. 
the  mean  effective  size  and  IWP.  are  incorporated 
m  radiation  calculations.  The  analysis  data  for 
winds,  temperature,  and  specific  humidity  at  12Z 
1  July  1979  are  ii.sed  as  initial  values  to  perform 
the  large  scale  prediction  for  cloud  fields. 


\Ve  'hall  confine  our  duciission  to  the  effects 
of  ihe  mean  e/'fective  si/e  of  ice  crvstals  on  large- 
'Cale  cloud  'imulations.  The  interactions  and  feed¬ 
backs  of  cloud  microphvsics,  radiative  transfer, 
and  the  cloud  model  are  shown  in  Fig.  4.  The 
model  predicts  IW'C  land  LWC).  cloud  cover,  tem¬ 
perature.  and  specific  humidity.  The  bulk  terminal 
velocity  and  mean  effective  size  are  defined  based 
on  the  cloud  microphvsics  parameterizations.  The 
IWC  together  with  the  model  thickness  give  the 
IWP.  The  terminal  velocity  is  directly  related  to 
IWC  by  means  of  parameterization  and  is  inter¬ 
active  in  the  cloud  model.  .-Xlthough  we  have  de¬ 
veloped  a  theoretical  relationship  between  IWC 
and  O,.  this  relationship  has  yet  to  be  verified.  In 
this  study,  we  shall  treat  D,  and  IWP  as  two  in¬ 
dependent  parameters  to  investigate  the  sensitivity 
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of  D,  on  the  Minulutions  of  large-scsilc  cloud  fields. 
Light  scattering  and  absorption  calculations  for 
cloud  particles  require  onU  information  on  par¬ 
ticle  si/e  distribution,  which  is  parameteri/ed  in 
terms  of  /),.  For  radiative  transfer  calculations, 
the  optical  depth  \alue  is  also  needed.  As  shown 
in  the  preceding  sections,  the  optical  depth  can  be 
e.xpressed  in  terms  of  D,  and  IWT.  Temperature 
and  humidity  are  inputs  in  radiative  transfer  cal¬ 
culations.  Lastly,  cloud  cover  from  the  model  is 
required  to  compute  the  average  radiative  heating 
in  a  model  grid  area.  Radiative  heating  affects 
temperature  directls  via  the  thermodsnamic  heat 
equation,  while  the  bulk  terminal  velocity  affects 
1\N'C  (and  LWC)  directly  through  the  governing 
equations  for  IWC.  LWC  and  precipitation. 

Comparisons  of  the  zonally  averaged  total 
cloud  cover  and  outgoing  longwave  radiation 
(OLR)  predicted  from  the  model  and  deri\ed  from 
obser\  ations  are  shown  in  Figs.  5a  and  5b.  The 
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Fip  5  C omparis.in'  ol  ihr  7on.ill>  averaged  (al  lolal  cloud 
cover  and  (hi  Ol.R  pieduted  from  ihc  model  and  derived 
from  .IDNFPH  and  ERH  for  Julv  |4''y  The  model  rc-ull'- 
are  for  averages  over  a  24  h  period 


predicted  civ’'ud  cewcis  arc  averages  of  the  re  .uUs 
computed  at  yf  lime  steps  over  a  24  hour  period. 
The  observed  cloud  covers  are  taken  from  Hen- 
derson-Sellers  ( 19K6)  based  on  the  3DNEPH  data 
base  for  July  1979.  It  is  well  known  that  the 
3DNEPH  under-  and  over-estimates  cloud  cover 
in  Arctic  and  Antarctic  regions,  respectively,  for 
the  summer  season  of  the  Northern  Hemisphere. 
During  this  period,  there  are  persistent  stratus 
clouds  over  the  Arctic  region.  However,  the 
threshold  niethod  utilized  in  3DNEPH  has  not 
been  able  to  distinguish  between  the  cold  surface 
and  clouds  aloft.  As  a  result,  because  of  cold  tem¬ 
peratures  in  the  .Antarctic  region  during  the  north¬ 
ern  summer,  the  threshold  technique  overesti¬ 
mates  cloud  cover.  The  model  underestimates 
cloud  cover  in  the  tropics  because  of  the  nature 
of  the  present  3-D  global  cloud  model,  w  hich  was 
developed  primarilv  for  stratiform  clouds.  It  ap¬ 
pears  that  improvements  could  be  made  if  a  cu¬ 
mulus  convection  program  were  added  in  the 
model.  Nevertheless,  there  are  general  similarities 
between  the  predicted  and  observed  cloud  covers. 
The  OLR  fields  computed  from  the  model  com¬ 
pare  reasonably  well  with  the  monthly  averaged 
data  from  ERB  for  July  1979.  Discrepancies  be¬ 
tween  the  tw  o  may  be  e.xplained  by  the  cloud  cover 
patterns  described  previously.  In  summary,  the 
similarities  between  the  model  results  and  obser¬ 
vations  in  terms  of  cloud  cover  and  OLR  are  en¬ 
couraging  in  view  of  the  simplicity  of  the  cloud 
model  structure.  Sensitivity  experiments  are  car¬ 
ried  out  to  investigate  the  relative  importance  of 
the  mean  effective  size  on  the  predicted  large-scale 
cloud  field. 

The  cloud  model  predicts  IWC  for  high  clouds 
and  LWCs  for  middle  and  low  clouds.  Mean  ef¬ 
fective  sizes  of  50  and  75  pm  are  used  in  large- 
scale  cloud  simulations.  The  value  of  50pm  cor¬ 
responds  to  the  mean  effective  size  for  a  typical 
cirrostratus.  Show  n  in  Figs.  6a,  6b  and  6c  are  the 
zonally  averaged  differences  in  the  IR,  solar  and 
net  heating  rates  predicted  by  the  model  for  a  24 
hour  period.  For  smaller  D^.  the  IR  emissivity  is 
greater  for  high  clouds,  leading  to  more  cooling 
at  high  cloud  tops.  However,  the  greater  high 
cloud  emissivity  produces  heating  at  high  cloud 
bases  and  increases  heating  rates  in  the  middle 
cloud  level  The  contour  lines  in  Fig.  6a  are  for 
each  0.3°C  dav  "  The  increase  of  high  cloud  top 
cooling  for  smaller  /),  is  on  the  order  of  about 
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0.3’C  day  ~  '  in  cloud  regions.  In  the  middle  cloud 
level,  heating  rates  increase  by  about  0.6°C  day  ' 
These  values  are  significant,  since  they  are  zonally 
averaged  quantities.  No  change  is  seen  in  the  heat¬ 
ing  rate  patterns  below  middle  clouds.  The  dif¬ 
ferences  in  solar  heating  rates  produced  by  the 
two  mean  effective  sizes  are  confined  in  the  North- 
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Fig.  6.  Zonallv  a\er.igcd  dilTerences  in  (a)  IR.  (b)  solar  and 
(c)  not  heating  rates  for  a  24  h  period  using  mean  effeetive 
sizes  lit  '0  .ind  '5|im  m  numeneal  experiments 


ern  Hemisphere  because  the  July  condition  is  used. 
.A  smaller  D^.  retlects  more  solar  radiation  and 
increases  solar  heating  in  high  cloud  top  regions. 
At  the  same  time,  it  reduces  the  solar  fiu.xes  avail¬ 
able  in  the  middle  cloud  level.  The  contour  lines 
in  Fig.  6b  are  for  each  0.02°C  day  which  is  15 
times  smaller  than  the  IR  heating  cooling  rates. 
Figure  6c  shows  the  net  heating  differences.  The 
patterns  are  similar  to  the  IR  heating, cooling  pat¬ 
terns.  The  e.xception  is  the  reduction  of  cooling 
above  high  cloud  top  regions  in  the  Northern 
Hemisphere. 

Figures  7a-7c  show  the  effects  of  the  mean  ef¬ 
fective  size  on  the  simulation  of  temperature,  cloud 
cover  and  IWC.  Because  of  more  cooling  pro¬ 
duced  by  a  smaller  D,.  above  the  high  cloud  top, 
temperatures  decrease  there  (Fig.  7a).  The  maxi¬ 
mum  decrease  is  on  the  order  of  about  0.2-0.4°C. 
However,  temperature  increases  in  the  middle 
cloud  level  by  about  0.4°C  across  the  latitudes.  In 
response  to  temperature  decreases  above  the  high 
cloud  top,  high  cloud  cover  increases.  The  increase 
of  high  cloud  cover  is  on  the  order  of  4%  across 
most  latitudes  (Fig.  7b).  The  use  of  different  mean 
effective  sizes  of  ice  crystals  significantly  affects 
the  calculation  of  middle  cloud  cover  due  to  tem¬ 
perature  increases.  E.xcept  in  the  vicinity  of  .30°S 
where  the  occurrence  of  clouds  is  minimal,  middle 
cloud  cover  is  reduced  by  3  to  4%  between  60°N 
and  60°S  latitudes  because  of  warmer  tempera¬ 
tures  produced  by  IR  heating.  In  terms  of  IWC 
for  high  clouds  (Fig.  7c).  differences  produced  by 
the  two  simulations  are  rather  small  because  IWC 
is  largely  controlled  by  the  bulk  terminal  velocity, 
which  is  interactive  in  the  model.  Variations  in  the 
mean  effective  size  do  not  affect  the  bulk  terminal 
velocity  in  the  present  experiments.  It  appears  that 
IWC  variations  are  primarily  governed  by  micro¬ 
physical  processes.  If  the  mean  effective  size  could 
be  made  interactive  through  IWC,  radiative  per¬ 
turbations  vvould  be  important  in  the  generation 
of  IWC. 

The  preceding  experiments  have  been  carried 
out  to  explore  the  importance  of  the  mean  effective 
size  of  ice  crystals  in  large-scale  cloud  simulations. 
Based  on  experimental  results,  small  mean  effec¬ 
tive  size,  through  radiative  heating  perturbations, 
would  increase  high  cloud  cover  and.  at  the  .same 
time,  reduce  middle  cloud  cover.  It  is  commonly 
recognized  that  the  prime  importance  of  high 
clouds  is  the  downward  emission  of  IR  radiation 
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Fig  7  Zonall)  a\crdgcj  dilTcrenccs  in  (a I  lenipcr.ilurc,  (b) 
cloud  cover  and  Ic)  cK'ud  IW  C’  for  a  24  h  period  using  mean 
effecliu'  sues  of  50  ansi  ’spm  in  numerical  cvperimcnts 

(greenhouse  effect)  For  middle  clouds,  it  is  the 
reflection  of  solar  radiation  (solar  albedo  effect) 
The  increase  of  high  cloud  cover  and  the  decrease 
of  middle  cloud  cover  would  lead  to  the  ampli¬ 
fication  of  the  greenhouse  effect  Small  mean  ef¬ 
fective  sizes  arc  normalK  considered  to  he  siL'nif- 


have demonstrated  in  this  studs  that  the  prevailing 
effects  of  small  mean  effective  sizes  are  in  the  pi  o- 
duction  of  IR  heating  and  cooling.  The  IR  healing 
and  cooling  in  the  atmosphere  produced  by  small 
mean  effective  sizes  are  much  more  significant 
than  the  solar  heating  counterpart. 


5.  Conclusions 

Parameterization  programs  for  cloud  microphys- 
ics  and  radiative  transfer  involving  ice  clouds  hav  e 
been  developed  in  conjunction  with  a  3-D  global 
cloud  model.  We  show  that  the  mean  effective  size 
and  IWP  are  two  variables  that  are  related  to  the 
optical  depth.  The  mean  effective  size  appears  suf¬ 
ficient  to  define  the  ice  crystal  size  distribution 
with  respect  to  light  scattering  and  absorption  cal¬ 
culations.  We  illustrate  the  importance  of  the 
mean  effective  size  in  radiative  transfer  calcula¬ 
tions  in  terms  of  solar  reflectivity.  IR  emissivity 
and  radiative  heating.  For  a  given  IW'P.  smaller 
mean  effective  sizes  reflect  more  solar  radiation, 
emit  more  IR  radiation  and  enhance  net  radiative 
heating  cooling  at  the  cloud  top  and  bottom  than 
larger  sizes. 

A  3-D  global  cloud  model  has  been  used  to 
investigate  the  effects  of  the  mean  effective  size  of 
ice  crystals  on  the  simulation  of  cloud  cover  and 
IW’C.  This  model  has  been  verified  through  com¬ 
parisons  of  the  predicted  cloud  cover  and  OLR 
with  monthly  averaged  values  derived  from 
3DNEPH  amd  ERB.  High  clouds  that  are  com¬ 
posed  of  smaller  mean  effective  sizes  would  pro¬ 
duce  more  cooling  at  the  cloud  top  and  more  heat¬ 
ing  at  the  cloud  base,  than  high  clouds  composed 
of  larger  sizes.  As  a  result  of  these  heating  cooling 
patterns,  more  high  cloud  cover  and  less  middle 
cloud  cover  are  simulated  from  the  model.  Using 
the  mean  effective  sizes  of  50  and  75  gm,  a  ~  4°o 
increase  and  3^%  decrease  are  shown  for  high 
and  middle  cloud  covers,  respectively.  The  radi¬ 
ative  heating  cooling  patterns,  however,  do  not 
affect  IWC  for  high  clouds  or  LW'C  for  middle 
clouds  because  the  mean  effective  size  is  prescribed 
in  the  model  If  a  relationship  could  be  established 
between  the  mean  effective  size  and  IWC,  radiative 
healing  through  the  mean  effective  size  would  in¬ 
fluence  the  prediction  of  IWC  as  well  as  LW'C  for 
middle  clouds 
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On  the  basis  of  the  preceding  diseu^si^'n,  it 
uould  be  important  to  develop  a  relationship  be¬ 
tween  the  mean  effeetive  size  of  iee  ersstals  and 
IWC  from  aircraft  experimental  data.  Second,  ra¬ 
diative  transfer  parameterizations  that  encompass 
the  mean  effective  size  developed  in  Section  .1 
should  be  verified  through  concurrent  radiation 
and  cloud  iidi^iophysics  experiments  involv  ing  ice 
clouds,  bdoally,  the  3-D  global  cloud  model  may 
be  improved  by  incorporating  time  varying  wind 
fields  to  study  their  effects  on  the  simulations  of 
cloud  parameters  in  the  context  of  cloud  micro¬ 
physics  radiative  transfer  parameterizations.  The 
program  may  be  accomplished  by  using  a  separate 
and  parallel  GC.M  lo  produce  the  required  wind 
fields. 
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ABSTRACT 


Maximum  infrared  polarization  signature  up  to  ~1%  is  predicted  in  tropical 
subvisual  cirrus  involving  randomly  oriented  ice  crystals  based  on  radiative 
transfer  calculations. 
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In  the  development  of  target  detection  based  on  infrared  signature  from 
space,  the  air  and  the  ground,  the  effects  of  cirrus  clouds,  particularly 
subvisual  cirrus,  on  infrared  transmission  must  be  understood.  The  radiative 
properties  of  cirrus  clouds  have  been  investigated  by  Takano  and  Liou^-^  using 
the  scattering  and  absorption  data  derived  from  hexagonal  ice  crystals.  In  these 
papers,  the  radiative  transfer  method  has  been  specifically  developed  for  solar 
wavelengths  in  which  thermal  emission  can  be  neglected.  Liou  et  al.^  have 
extended  the  radiative  transfer  methodology  to  the  thermal  infrared  wavelengths 
in  order  to  study  the  transmission  properties  of  thin  cirrus  clouds.  In  Takano 
and  Liou^,  we  have  demonstrated  that  there  are  important  polarization  signatures 
of  reflected  sunlight  from  clouds.  These  signatures  could  be  utilized  to 
differentiate  between  spherical  and  nonspherical  particles,  as  well  as  to 
distinguish  among  various  nonspherical  particles. 

However,  the  information  concent  of  the  polarization  from  cirrus  clouds  in 
the  thermal  infrared  wavelengths  has  not  been  explored.  The  contribution  to  sky 
polarization  from  hexagonal  and  oriented  ice  crystals  could  be  significant.  It 
is  the  objective  of  this  note  to  present  some  results  regarding  the  radiance  and 
polarization  pattern.^  for  the  10  /im  wavelength  in  a  number  of  cirrus  cloudy 
atmospheres . 

The  single-scattering  properties  for  randomly  oriented  hexagonal  ice 
crystals  in  the  infrared  wavelengths  have  been  computed  using  the  program 
developed  by  Takano  and  Liou'.  Figure  la  shows  the  phase  function  for  hexagonal 
ice  crystals  with  an  aspect  ratio  (L/2a)  of  120/60  (/im/^m)  at  the  10  fxm 
wavelength,  where  L  is  the  length  and  2a  is  the  width  of  a  hexagonal  crystal. 
The  complex  refractive  index  of  ice*  used  in  this  computation  is  1 . 1991 - iO . 051 . 
For  comparison,  the  phase  function  for  the  area-equivalent  Mie  spheres  is  also 
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shovm  in  this  figure.  The  peaks  on  the  ice  crystal  phase  function  at  the 
scattering  angles  of  14'  and  27'  correspond,  respectively,  to  inner  and  outer 
haloes,  which  are  at  different  angles  than  the  standard  22'  and  46'  haloes  in  the 
visible  because  of  the  differences  in  refractive  indices.  Ice  crystals  scatter 
more  light  in  backward  directions  than  spheres.  Figure  lb  shows  the  degree  of 
linear  polarization  due  to  single  scattering  of  hexagonal  crystals  and  spheres. 
Generally,  the  degree  of  linear  polarization  lor  hexagonal  crystals  is  smaller 
than  that  for  ice  spheres.  Using  the  single-scattering  and  polarization  data, 
computations  of  the  transfer  of  polarized  infrared  radiation  in  cirrus  cloudy 
atmospheres  have  been  carried  out.  Under  plane  parallel  conditions,  radiance  and 
polarization  patterns  are  azimuthally  independent.  They  are  only  functions  of 
the  optical  depth  and  zenith  angle.  In  the  computations,  the  surface  is  assumed 
to  emit  unpolarized  radiation  isotropically. 

We  have  used  tropical,  mid-latitude  winter,  and  arctic  winter  atmospheric 
profiles  in  performing  the  radiance  and  polarization  calculations.  The  cloud  and 
surface  temperatures  for  these  profiles,  which  are  of  prime  importance  in 
infrared  radiative  transfer  calculations,  are  listed  in  Table  1.  The  heights  of 
cirrus  clouds  are  also  included  in  this  table.  We  used  a  columnar  ice  crystal 
shape  with  an  aspect  ratio  (length/diameter)  of  120/60  (pm/pm)  whose  single 
scattering  properties  at  10  /an  are  shown  in  Fig.  1. 

Figure  2a -d  shows  the  upward  radiances  at  the  cloud  top  and  the  downward 
radiances  at  the  cloud  base.  Limb  darkening  and  limb  brightening  are  evident  in 
upward  and  downward  patterns,  respectively.  These  patterns  can  also  be  seen  in 
the  window  radiances^.  When  the  cloud  optical  depth  is  small  (Tj.-1),  the  upward 
radiance  for  zenith  angles  less  than  -60®  are  close  to  those  emitted  from  the 
surface.  For  zenith  angles  larger  than  -60',  the  upward  radiance  is  reduced  due 
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to  cloud  absorption  of  the  radiation  emitted  from  the  surface.  However,  the 
downward  radiance  is  enhanced  due  to  the  combined  effects  of  cloud  emission  and 
the  cloud  reflection  of  surface  radiation.  The  upward  radiances  decrease  with 
the  increase  of  optical  depth  in  any  models,  because  cirrus  clouds  absorb  the 
radiation  emitted  by  the  surface.  In  contract,  the  downward  radiances  increase 
with  the  increase  of  optical  depth  in  any  models  because  of  the  emission  of 
cirrus.  When  the  cloud  optical  depth  is  large  (r^-A),  the  downward  radiances  for 
zenith  angles  less  than  -60®  are  close  to  those  emitted  from  the  cloud.  As  shown 
in  this  figure,  radiance  patterns  depend  on  the  atmospheric  temperature  profile 
used  in  the  calculation.  For  r^-l,  the  upward  radiances  increase  according  to  the 
surface  temperature.  Exceptions  are  for  the  tropical  case  involving  zenith 
angles  greater  than  about  70®  because  of  the  lower  temperature  of  the  tropical 
anvil.  For  optically  thick  clouds  with  the  upward  radiances  are  largely 
produced  by  the  cloud-top  emission.  The  downward  radiances  depend  primarily  on 
the  cloud  temperature. 

The  degree  of  linear  polarization®  patterns  corresponding  to  Fig.  2a-d  is 
illustrated  in  Fig.  3a-d,  Because  of  the  differences  in  upward  and  downward 
radiances,  the  linear  polarization  patterns  are  also  different  at  the  cloud  top 
and  base.  It  is  noted  that  there  is  a  significant  difference  in  linear 
polarization  patterns  between  thermal  infrared  and  solar  wavelengths.  The 
magnitude  of  linear  polarization  displayed  ^n  Fig.  3,  that  is  generally  less  than 
1%,  is  much  smaller  than  that  in  the  case  of  solar  radiation  (see  Fig.  8  in 
Takano  and  Liou^) .  This  is  in  view  of  the  fact  that  the  surface  and  a  small 
volume  of  cloud  are  assumed  to  emit  unpolarized  thermal  radiation  isotropically. 
We  shall  consider  the  optically  thin  case  (r^-l)  in  Fig.  3a-b  in  which  the  effect 
of  single  scattering  events  on  polarization  Is  most  pronounced.  In  this  case. 
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radiation  emerged  from  the  cloud  top  is  dominated  by  the  emission  from  the 
surface  and  largely  unpolarized  for  zenith  angles  less  than  60®.  For  larger 
zenith  angles,  positive  polarization  may  be  produced  by  single  scattering  of  the 
emitted  unpolarized  radiation  from  the  surface  as  illustrated  in  Fig.  lb.  At  the 
cloud  base,  positive  linear  polarization  is  produced  also  due  to  the  single 
scattering  of  the  emitted  unpolarized  radiation  from  the  surface.  When  is 
1/8,  downward  polarization  can  reach  more  than  1%  at  ^“60®.  When  the  cloud 
optical  depth  increases,  polarization  that  is  produced  by  single  scattering  is 
generally  reduced.  However,  total  polarization  does  not  monotonically  depend  on 
optical  depth,  because  of  the  combined  effect  of  radiances  from  the  cloud  and  the 
surface . 

From  Fig.  3a  and  c,  it  can  be  generally  concluded  that  the  upward 
polarization  at  the  top  of  clouds  composed  of  randomly  oriented  ice  crystals  is 
small,  and  for  all  practical  purposes,  the  polarization  can  be  neglected.  In  the 
case  of  downward  components,  polarization  emerged  from  the  tropical  cirrus  is 
greater  than  that  emerged  from  the  midlatitude  and  arctic  cirrus  by  about  one 
order  of  magnitude  as  shown  in  Fig.  3b  and  d.  This  feature  can  be  explained  as 
follows:  First,  consider  only  a  cirrus  cloud  without  an  underlying  surface.  The 
degree  of  linear  polarization  at  cloud  boundaries  is  slightly  negative  and 
independent  of  the  cloud  temperature.  This  is  because  the  effect  of  cloud 
temperature  is  canceled  out  in  the  denominator  and  numerator  of  -Q/1 .  Second, 
we  add  a  surface  below  a  thin  cirrus.  The  total  polarization  is  largely  affected 
by  the  single  scattering  of  surface -emitted  radiation  by  ice  crystals. 
Generally,  single-scattered  radiation  has  positive  linear  polarization  as  pointed 
out  previously.  The  temperature  difference  between  the  surface  and  the  cloud  in 
the  case  of  tropical  cirrus  is  greater  than  that  for  the  midlatitude  and  arctic 
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cirrus.  As  a  result,  the  tropical  cirrus  shows  more  Intense  positive 
polarization  than  the  midlatitude  and  arctic  cirrus.  Downward  polarization  for 
water  clouds  will  also  be  small  due  to  smaller  differences  between  surface  and 
cloud  temperatures. 

Radiance  and  polarization  patterns  have  also  been  calculated  for  Ice  plates 
with  an  aspect  ratio  of  20/20  (pm/pm).  Computed  results  are  close  to  those  for 
ice  columns  with  an  aspect  ratio  of  120/60  (pm/pm).  The  radiance  and  degree  of 
linear  polarization  at  the  cirrus  cloud  boundaries  depend  more  significantly  on 
the  difference  between  the  surface  and  cloud  temperatures  than  on  the  cloud 
particle  size  and  shape.  Downward  infrared  polarization  at  the  cloud  base  in  the 
case  of  tropical  cirrus  can  reach  more  than  1%,  which  is  greater  than  that  for 
the  cases  of  midlatitude  and  arctic  cirrus  by  about  one  order  of  magnitude.  At 
t-he  wavelength  of  10  pm,  optical  depths  of  air  molecules  and  aerosols  are  small. 
Thus,  detection  of  the  tropical  subvisual  cirrus  may  be  possible  by  measuring 
downward  polarization  at  zenith  angles  <f  60* -70*  from  the  surface.  This 
possibility  will  be  explored  in  our  future  research  efforts.  The  effects  of 
lorizontal  orientation  of  ice  crystals  and  the  associated  anisotropic  emission 
on  infrared  polarization  are  also  subjects  for  further  investigation. 
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Figure  Captions 


Fig.  1 


Fig.  2 


Fig.  3 


(a)  Scattering  phase  function  and  (b)  degree  of  linear  polarization 
for  randomly  oriented  hexagonal  ice  crystals  of  L/2a  -  120/60 
itim/tim)  and  the  area  equivalent  Mie  spheres  at  A  -  10  pm. 

Radiances  as  a  function  of  zenith/nadir  angles  on  the  boundaries  of 
cirrus  clouds  in  the  three  atmospheric  profiles  in  the  cases  of  (a) 
and  (b)  r^-l  and  (c)  and  (d)  Tj-A. 

Degree  of  linear  polarization  as  a  function  of  zenith/nadir  angles 
on  the  boundaries  of  cirrus  clouds  in  the  three  atmospheric  profiles 
in  the  cases  of  (a)  and  (b)  r^-l  and  (c)  and  (d)  Tc-4. 


Table  1.  Cloud  and  surface  temperature  and  cloud  heights  for  three 
atmospheric  profiles. 


Profile 

Cloud 

Height  (km) 

Cloud 

Temperature ,  Tj.  (K) 

Surface 

Temperature,  T,  (K) 

Tropical 

-15.0 

200 

300 

Midlatitude  Winter 

-8.0 

230 

272 

Arctic  Winter 

-6.5 

227 

257 
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